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The  charge  transport  process  in  InGaAs/InAlAs 
multiquantum-well  (MQW)  heterostructure  p-i-n  photodiodes 
is  studied  experimentally  and  theoretically.  Current- 
voltage  (I— V) , capacitance-voltage  (C-V) , S-parameter,  and 
noise  measurements  are  performed  at  room  temperature  over 
a wide  range  of  reverse  bias.  The  devices  under  study 
consist  of  alternating,  undoped,  lattice  matched  InGaAs 
and  InAlAs  layers  sandwiched  between  a heavily  doped  n and 
p layer  of  InGaAs.  To  determine  the  dependence  of  charge 
transport  phenomena  on  quantum-well  width,  a set  of  three 
different  quantum-well  layers  of,  respectively,  30,  90, 
and  500  Angstrom  wide  were  used. 

Upon  illumination  by  a light-emitting-diode  (1350 
nm) , measurements  of  the  noise  spectra  were  performed  at 
frequencies  ranging  from  100  Hz  to  1 GHz.  In  the  frequency 
range  of  200  MHz  and  1 GHz,  the  90-90  and  500-500  Angstrom 
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devices  which  have  nontransparent  barriers  between 
neighboring  quantum -wells,  show  distinctive  sub-shot  noise 
level  at  low  bias.  This  sub-shot  noise  phenomenon  is 
interpreted  in  terms  of  the  trapping  and  detrapping  of 
electrons  in  undepleted  quantum-wells.  The  magnitude  of 
the  noise  is  proportional  to  1/N,  where  N is  the  number  of 
undepleted  quantum-wells.  The  dependence  of  noise  levels 
on  bias  indicates  that  quantum-wells  located  in  the  i- 
region  are  successively  and  gradually  depleted  with 
increasing  reverse  bias.  C-V  measurements  support  this 
finding. 

To  explain  the  sub-shot  noise  levels,  a MQW  trapping 
noise  model  was  developed  employing  the  Transfer  Impedance 
method. 

The  excess  noise  observed  at  frequencies  below  200 
MHz  is  caused  by  the  trapping  and  detrapping  of  electrons 
in  hetero junction  interface  states. 

Current-voltage  measurements  performed  under 
illumination  indicated  that  the  quantum  efficiency  of  the 
photodiode  is  bias  dependent.  We  established  that  the 
photocurrent  is  controlled  by  field  assisted  injection  of 
optically  generated  electrons  from  the  low  bandgap 
absorption  region  into  the  intrinsic  region. 

Finally,  nonparabolicity  of  the  conduction  band  was 
studied  and  a new  formula  for  the  ionization  threshold 
energy  was  derived  and  compared  with  experimental  values. 
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CHAPTER  I 
INTRODUCTION 

In  recent  years  there  have  been  a number  of 
remarkable  advances  in  the  fabrication  of  crystalline 
semiconductor  heterostructures  made  possible  by  molecular 
beam  epitaxy  (MBE)  and  metal-organic  chemical  deposition 
(MOCVD) . These  epitaxial  growth  techniques  allow  for  the 
precise  control  of  both  the  thickness  of  heterostructure 
layers  and  the  composition  of  compound  materials  and  make 
bandgap  engineering  possible.  A new  generation  of  devices 
with  artificial  band  profiles  and  structures  is  emerging 
from  this  approach. 

Among  novel  structures,  multiquantum-well  (MQW) 
hetero junctions  and  superlattices  are  of  great  interest 
because  they  are  the  simplest  examples  of  man-made 
structures  which  allow  for  the  modification  of  electronic 
and  optical  properties.  Applications  as  well  as 
fundamental  studies  on  these  structures  are  being  worked 
on  by  researchers  in  many  different  areas. 

In  the  area  of  fiber-optic  communications, 
considerable  efforts  have  been  devoted  to  the  development 
of  low-noise,  high-gain  avalanche  photodiodes  (APDs) . As 
suggested  by  McIntyre,  excess  noise  in  APDs  is  small  when 
the  electron-to-hole  impact  ionization  coefficient  ratio, 
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a/p,  is  very  large  or  very  small.  However,  most  III-V 
semiconductors  of  interest  have  nearly  equal  electron  and 
hole  ionization  coefficients.  Therefore,  there  is  a strong 
desire  to  artificially  enhance  this  ratio. 

Chin  et  al . [1]  first  proposed  a method  for 
artificially  enhancing  the  impact  ionization  coefficient 
ratio  by  employing  a multilayered  GaAs/AlGaAs  structure. 
Capasso  et  al . [2]  and  Juang  et  al . [3]  reported  the 
experimental  results  of  an  enhanced  a/p  ratio  in 
GaAs/AlGaAs  multiquantum-well  structures. 

For  fiber-optic  communications,  wavelengths  of  1.3  pn 
or  1.55  |im  should  be  used  to  minimize  the  attenuation  of 
light  transmission  and  dispersion  in  optical-fibers.  Among 
the  materials  sensitive  to  these  wavelength,  the 
Ino . 53Gao . 47AS  ternary  compound  in  conjunction  with 
Irio.52Alo.48As  is  most  widely  used  to  make  hetero  junctions . 
Capasso  et  al.  [4]  experimentally  investigated 
Ino.53Gao.47As/Ino.52Alo.48As  multiquantum-well  structures  and 
reported  the  enhanced  a/p  ratio  at  low  temperature. 

One  of  our  initial  research  objectives  was  to  measure 
the  a/p  ratios  from  noise  measurements  in  several  MQW 
structures  with  different  dimensions.  However,  it  turns 
out  that  in  this  material  and  in  MQW  structures,  Zener 
tunneling  current  easily  dominates  with  the  result  that  we 
could  never  obtain  a high  enough  field  to  have  band-to- 
band  impact  ionization.  Allam  et  al . [5]  also  reported 
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that  strong  resonant  Zener  tunnelling  occurs  in  these 
structures  at  room  temperature.  Without  finding  a new  way 
to  reduce  the  Zener  tunnelling  current  an  InGaAs/InAlAs 
MQW  structure  is  not  suited  for  APD  applications.  This 
disappointing  situation  was  not  recognized  until  the 
experiments  were  actually  performed.  Therefore,  a more 
general  treatment  in  terms  of  charge  transport  is  needed 
to  investigate  this  novel  MQW  structure  before  considering 
its  applications. 

The  charge  transport  in  our  structures  takes  place  in 
the  direction  perpendicular  to  the  hetero junction 
interfaces.  In  semiconductor  transport  problems  the 
following  four  situations  may  arise.  Note  that  a 
combination  of  more  than  one  situation  may  exist. 

1)  The  carrier  density  increases  due  to 
multiplication . 

2)  The  carrier  density  decreases  due  to 
recombination . 

3)  The  carrier  flow  is  slowed  down  due  to  trapping  or 
scattering. 

4)  The  carrier  flow  speeds  up  due  to  higher  mobility 
and  less  scattering. 

In  our  structure  the  first  situation  is  excluded  as 
explained  above.  The  second  case  does  not  apply  to  the 
depleted  section  of  quantum-well  layers.  The  third  case  is 
most  probable  to  occur  in  the  MQW  structures  we 
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investigated.  Finally,  the  fourth  situation  usually 
applies  to  the  situation  where  carrier  transport  is  in  a 
direction  parallel  to  the  heterostructure  interfaces,  as 
it  occurs  for  example  in  modulation  doped  FETs. 

The  main  objective  of  this  dissertation  is  to  study 
the  charge  transport  noise  in  MQW  heterostructure  p-i-n 
diodes.  The  noise  measurement  tool  was  extensively 
employed  to  investigate  the  microscopic  charge  transport 
mechanisms  of  the  p-i-n  photodiodes.  To  measure  the  high 
frequency  noise,  a new  measuring  set-up  was  built.  The 
noise  spectrum  was  measured  over  a wide  frequency  range 
from  100  Hz  to  1 GHz.  In  addition  to  noise  measurements, 
conventional,  basic  measurements  such  as  current-voltage 
( I— V) , capacitance-voltage  (C-V)  measurements  and 
impedance  (S-parameter ) measurements  were  performed  to 
characterize  the  devices  under  study.  As  an  unexpected 
finding,  we  observed  distinctive  sub-shot  noise  levels  at 
low  bias  which  to  the  author's  knowledge  have  never  been 
observed  before.  In  Chapter  VI,  this  sub-shot  noise 
phenomenon  is  identified  and  interpreted  in  terms  of 
trapping  and  detrapping  in  the  quantum-wells.  To  explain 
the  results  of  our  noise  measurements  we  developed  a 
simplified  trapping  and  detrapping  model  and  a qualitative 
theory. 

The  organization  of  this  thesis  is  as  follows. 
Chapter  II  presents  a review  of  conventional  avalanche 
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photodiodes.  The  details  of  the  impact-ionization 
mechanism  and  device  noise  are  presented.  In  addition, 
various  design  considerations  of  an  APD  are  discussed  as 
they  relate  to  low-noise  performance. 

In  Chapter  III,  the  expression  for  the  excess  noise 
factor  of  an  ideal  n-stage  staircase  multiplier  is 
presented  with  some  proofs. 

Chapter  IV  discusses  nonparabolicity  as  a tool  in 
band-gap  engineering.  A new  formula  for  the  threshold 
energy  of  electron-electron-hole  impact  ionization  for  the 
case  of  a parabolic  valence  band  and  a simple  Kane  type, 
nonparabolic  conduction  band  is  presented. 

In  Chapter  V,  the  theory  for  carrier  trapping  and 
detrapping  in  multiquantum-wells  and  at  hetero junction 
interfaces  is  developed.  As  a result,  an  expression  for 
the  sub-shot  noise  factor  is  found. 

Chapter  VI,  the  highlight  of  this  dissertation, 
describes  the  multiquantum-well  device  structures  and 
various  measurement  schemes  used  in  this  study. 
Particularly,  the  noise  measurements  are  discussed  in 
detail.  Experimental  results  are  presented  with 
qualitative  discussions. 

Finally,  conclusions  are  presented,  followed  by  a 
summary  of  this  work  and  some  recommendations  for  future 


research. 


CHAPTER  II 

REVIEW  OF  CONVENTIONAL  AVALANCHE  PHOTODIODE  OPERATION 

2.1.  Impact  Ionization 

2.1.1.  Definition  of  the  Impact  Ionization  Coefficient 

When  the  electric  field  is  constant  in  the  avalanche 
region  and  a spatial  steady  state  is  achieved,  the  impact 
ionization  coefficient  is  defined  as  the  probability  per 
unit  length  along  the  direction  of  the  field  that  a 
carrier  impact  ionizes.  Equivalently,  this  definition 
reduces  to  the  reciprocal  of  the  average  distance  <x> 
along  the  direction  of  the  field  between  two  consecutive 
ionizing  collisions. 

The  ionization  coefficients  are  calculated  as  follows 

[6] 


1 

lvn|tn 


1 

lvplxp 


(2.1) 


where  vn  and  vp  are  the  electron  and  hole  drift 
velocities,  respectively,  and  xn  and  xp  are  the  times  for 
an  electron  and  hole  to  create  a pair.  Hence,  the 
ionization  rates  are  identified  with  the  mean  free  paths 
between  ionizing  collisions.  The  above  definition  is  not 
valid  for  situations  of  non-constant  electric  field  or  for 


6 


7 


an  ionization  coefficient  which  depends  on  position. 
Thornber  [7]  has  introduced  the  following  more  general 
definition  in  terms  of  the  number  of  electron-hole  pairs 
created  per  unit  time  per  unit  volume  normalized  to  the 
incident  carrier  density; 


where  £ denotes  the  electric  field,  Eth  denotes  the 
threshold  energy,  and  P (Eth)  represents  a probability 
function  which  depends  on  threshold  energy  and  electric 
field.  Since  P is  a strongly  increasing  function  of  the 
electric  field,  as  we  will  show  in  Eq.  (2.9),  P approaches 
unity  at  high  fields.  Then  the  average  distance  between 
consecutive  ionizing  collisions  becomes 


This  distance  is  the  minimum  distance  required  for  a 
carrier  to  create  an  electron-hole  pair. 

At  low  electric  fields,  where  P (Eth)  is  less  than 
unity,  this  distance  must  be  much  larger  than  the  minimum 
distance  and  is  expressed  as 


(2.2) 


a 1 ~ <x>  = (high  field) 


(2.3) 


Eth 


<x>  » 


e£ 


(low  field) 


(2.4) 
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2.1.2.  impact  Ionization  Mechanism 

Impact  ionization  is  an  essential  mechanism  in  the 
operation  of  avalanche  photodiodes  (APD)  and  impact- 
avalanche  transit  time  (IMP ATT)  diodes.  At  present  this 
mechanism  is  not  completely  understood  and  several 
competing  theories  are  available  [8,9,10]. 

Wolff  was  the  first  to  calculate  the  ionization  rate 
in  semiconductors.  He  applied  a method  borrowed  from  the 
gas  discharge  theory  to  solve  the  Boltzmann  transport 
equation.  He  assumed  that  electrons  experience  only  two 
types  of  collisions:  collisions  with  optical  phonons  and 
collisions  with  valence  electrons;  the  latter  is  producing 
impact  ionization. 

At  high  electric  fields,  he  assumed  that  the  energy 
gained  between  impact  ionizations  is  much  larger  than  the 
energy  lost  per  phonon  collision,  i.e., 

e'EX  » Ep  (2.5) 

where  Ep  is  the  energy  lost  per  phonon  collision  and  X is 
the  mean  free  path  between  impact  ionizations.  In  this 
situation  an  electron  can  undergo  many  phonon  collisions 
while  gaining  energy  from  the  electric  field.  Since  the 
electron  changes  its  direction  at  each  collision,  many 
phonon  collisions  serve  to  keep  the  velocity  distribution 
nearly  isotropic.  Hence,  using  only  the  first  two  terms 
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of  a Legendre  polynomial  expansion,  Wolff  approximated  the 
velocity  distribution  function  by  [8] 

F(v,e)  = fQ  (v)  + fi(v)cos0  (2.6) 

where  fQ (v)  represents  the  spherically  symmetric  part  of 
the  distribution  and  fi (v)  is  the  asymmetric  part  arising 
from  the  perturbations  caused  by  the  applied  electric 
field.  Since  9 is  the  angle  between  the  velocity  of  the 
carrier  and  the  electric  field,  the  cosG  term  represents 
the  net  contribution  to  the  distribution  function  in  the 
direction  parallel  to  the  electric  field. 

Using  this  approximate  distribution  function  and  a 
Maxwellian  approximation,  Wolff  solved  the  Boltzmann 
transport  equation  to  calculate  the  ionization  rate  with 
the  result 

a(£)  ~ exp  (-  ~)  (2.7) 

where  A is  a constant. 

On  the  other  hand  we  may  look  at  the  situation 

e£? i < Ep  (2.8) 

i.e.,  the  energy  gained  between  impact  ionization  is  less 
than  the  energy  lost  per  phonon  collision.  In  this 
situation  carriers  will  not  be  able  to  gain  the  high 
energies  required  for  impact  ionization  unless  they  escape 
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phonon  collisions.  These  carriers  are  the  "lucky"  or 
ballistic  carriers  first  considered  by  Shockley  [9] . The 
probability  that  a lucky  electron,  starting  at  zero 
energy,  ballist ically  gains  the  threshold  energy  for 
ionization  is  equal  to  the  probability  that  it  travels  for 
the  ionization  distance,  Eth/eE,  without  undergoing  a 
phonon  collision;  i.e.. 


P(Eth)  = exp  (-  ^k), 

elA 

from  which  Shockley  derived 

e£  E-t-h 

a = r; — exp  ( — ) 

Eth  p e£X 


(2.9) 


(2.10) 


Next,  Baraff  [10]  argued  that  the  distribution  function  of 
the  electron  energy  was  neither  a spike  in  the  direction 
of  the  field  (Shockley's  ballistic  theory)  nor  nearly 
isotropic  (Wolff) . He  said  that  these  two  aspects  were 
complementary  and  assumed  a velocity  distribution  of  the 
form 


f (v , COS0)  - A (v)  + B (v)  5 (1  - COS0)  (2.11) 

which  is  a mixture  of  a spherical  part  and  a "spike"  in 
the  direction  of  the  electric  field. 

By  using  this  distribution  function,  Baraff  solved 
the  Boltzmann  equation  and  showed  that  his  result 
contained  Shockley's  result  as  a low-field  limit  and 
Wolff's  result  as  a high-field  limit. 
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His  theory  is  fairly  successful  in  describing  the 
ionization  rates  in  Si  and  Ge,  but  does  not  explain  very 
well  the  ionization  rates  in  materials  where  strong  band 
structure  effects  exists,  as  in  III-V  compounds.  Also,  his 
theory  somewhat  obscures  the  detailed  physical  mechanism 
by  which  carriers  reach  the  ionization  energy. 

Shichijo  and  Hess  [11]  performed  a Monte  Carlo 
simulation  of  the  impact  ionization  phenomenon  in  GaAs  and 
provided  a very  clear  physical  picture  of  the  impact 
ionization  process.  They  explained  it  as  follows. 

Most  of  the  time,  an  electron  stays  around  the 
average  energy  of  0.7  ~ 1.2  eV  under  the  heating  influence 
of  an  electric  field  of  the  order  of  5 x 105  V/cm  and  the 
cooling  influence  of  phonon  collisions.  Occasionally, 
however,  one  of  these  electrons  will  escape  phonon 
collisions  and  move  up  to  higher  energies.  When  an 
electron  reaches  2.0  eV,  as  in  the  case  of  GaAs,  it  causes 
an  impact  ionization.  We  can  think  of  this  electron  as  one 
of  the  lucky  electrons  in  Shockley's  theory  [9],  and  of 
the  electrons  around  the  average  energy  as  the  diffusing 
electrons  in  Wolff's  theory  [8]. 

It  is  also  important  to  note  the  difference  between 
this  theory  and  Shockley's  theory.  Shockley's  "lucky" 
electrons  start  from  zero  energy,  escape  the  phonon 
scattering  completely,  and  impact  ionize.  The  ionizing 
electrons  in  the  Hess  theory  start  at  an  average  energy 
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and  reach  ionization  threshold  after  a few  scattering 
event  s . 

Since  this  theory  accounts  for  realistic  scattering 
mechanisms  and  band  structures,  it  can  be  applied  to  any 
semiconductor.  Shichijo  and  Hess  [11]  also  calculated  the 
electron-ionization  rate  in  GaAs  versus  crystal 
orientation  and  concluded  that  within  the  statistical 
uncertainty  of  their  calculations  (±  20  %) , no  orientation 
dependence  of  impact  ionization  for  electric  fields  larger 
than  3.3  x 106  V/cm  was  found.  They  explained  this  as 
follows:  At  high  electric  fields  and  at  steady  state, 
electrons  are  scattered  all  over  the  Brillouin  zone  before 
reaching  threshold,  so  any  "memory"  of  the  orientation  of 
the  electric  field  or  the  effect  of  the  anisotropy  of  the 
conduction  band  are  lost . 

Recently,  Ridley  [12]  has  developed  a simple 
analytical  theory  of  impact  ionization  (the  "lucky  drift 
model") , giving  results  that  are  in  excellent  agreement 
with  Baraff's  theory  [10].  The  basis  of  his  approach  is 
the  distribution  between  the  rates  of  momentum  relaxation 
Tm  and  energy  relaxation  xE.  If  we  have,  as  generally  is 
the  case,  Tm(x)  « te(E),  then  it  is  possible  for  electrons 
to  drift  in  an  electric  field  with  a drift  velocity 
limited  by  momentum-relaxation  collisions,  without  there 
being  in  the  same  period  any  significant  energy 
relaxation.  This  is  called  lucky  drift,  and  it  is 
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intermediate  between  Shockley's  ballistic  state  and 
Wolff's  equilibrium  state,  each  of  which  is  valid  in  a 
distinctive  time  domain  and  which  can  be  characterized  as 
follows  [12] : 

0 < t < xm(x);  ballistic  motion  regime 
Tm(T)  < t < te(E);  lucky  drift  regime 
te(E)  < t;  steady  state  regime 

2.1.3.  Impact  Ionization  Threshold  Energy 

Impact  ionization  can  occur  in  a semiconductor  when 
free  carriers  (electrons  and  holes)  are  accelerated  by  a 
high  electric  field  until  they  gain  sufficient  energy  to 
create  an  electron-hole  pair  by  the  generation  of  an 
electron  from  the  valence  band  to  the  conduction  band. 

The  minimum  energy  needed  for  impact  ionization  is 
the  ionization  threshold  energy  Eth*  The  value  of  the 
threshold  energy  depends  on  the  band  structure  of  the 
semiconductor.  The  conditions  for  an  impact  ionization 
process  to  occur  at  threshold  were  first  established  by 
Keldysh  [13].  Anderson  and  Crowell  [14]  developed  an 
algorithm  to  determine  the  threshold  condition  for  a 
general  type  of  band  structure.  Using  the  same  methods, 
Pearsall  [15]  obtained  the  relation  between  the  impact 
ionization  threshold  energy,  the  energy  band  gap,  and  the 
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effective  masses  of  the  carriers  (electrons  and  holes)  for 
GaAs  and  GaSb. 

In  the  following  we  will  review  the  above-mentioned 
papers  by  Anderson  and  Crowell  and  by  Pearsall  et  al . 
[14,15] . 

2. 1.3.1.  Conditions  for  Threshold. 

Fig.  2.1  shows  a general  band  structure  of  a direct 
bandgap  semiconductor  with  an  impact  ionization  process 
indicated.  An  initial  electron  in  the  conduction  band 
makes  a transition  from  state  i(Ei,ki)  to  state  l(Ei,ki), 
promoting  an  electron  to  state  2(E2,k2)  and  resulting  in  a 
hole  in  the  valence  band  state  3(E3,k3).  The  two  key 
requirements  in  order  for  impact  ionization  to  occur  are 
in  terms  of  electronic  quantities: 

(1)  conservation  of  energy,  i.e. 

Ei  + Ev (k3)  = Ec(ki)  + Ec (k2)  (2.12) 

(2)  Conservation  of  momentum 

hki  + hk3  = hki  + tik2  (2.13) 

As  Anderson  and  Crowell  [14]  have  shown,  the  threshold 
energy  is  determined  by  minimizing  the  energy  of  the 
initial  particle  with  respect  to  small  changes  in  the  wave 
vectors  of  the  particles  (carriers)  involved.  Satisfying 
the  restrictions  of  conservation  of  energy  and 
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Figure  2.1  Simple  three  band  E-k  diagram  showing 

the  electron-initiated  impact  ionization. 
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conservation  of  momentum,  we  find  for  the  threshold  energy 
Ei 

dEi  = 0 = dkrVkE(ki)  + dk2-VkE(k2)  - dk3-VkE(k3)  (2.14) 
and  for  a given  ki 

dki  = 0 = dki  + dk2  - dk3  (2.15) 

Substituting  Eq.  (2.15)  into(2.14)  and  recognizing 

VkE(k)/h  as  a group  velocity,  we  find 

dk2 • (v2  - vi)  + dk3- (vi  - v3)  = 0 (2.16) 

To  satisfy  Eq.  (2.16)  for  arbitrary  dk2  and  dk3,  one 
requires 

vi  = v2  = v3  (2.17) 

That  is,  a necessary  condition  for  an  impact-ionizing 
carrier  to  have  minimum  energy  is  that  all  final  carriers 
have  equal  group  velocities  or  lie  on  the  same  axis  in  E-k 
diagram  [14].  Note  that  this  axis  should  include  the 
initial  state. 

2. 1.3. 2.  Calculation  of  Ionization  Threshold  for  Simple 
Models. 

To  get  a functional  expression  showing  how  the 
threshold  energies  depend  on  the  parameters  of  the  band 
structure,  we  study  the  following  simple  models. 
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(a)  Two  Parabolic  Bands 

As  Fig.  2.1  shows,  we  consider  one  conduction  band 
with  effective  mass  me  and  valence  band  with  effective 
mass  mh;  me,mh  > 0.  Using  the  threshold  condition  Eq. 
(2.17),  we  obtain  the  following  relationship  between  the 
wave  vectors  of  the  final  particles  using  that  the  mass  of 
an  electron  in  the  valence  band  is  equal  to  - mh 

ki  = k2  = -k3^  (2.18) 

Also,  because  of  momentum  conservation, 

ki  = 2ki  - k3  = ki(2  + ^)  (2.19) 

me 

If  we  choose  the  bottom  of  the  conduction  band  as  our 
reference,  it  should  hold  that 


Ei  = 


,2.  2 

h kj 

2me 


fr2kj  mh  2 

2 me  me  Eth,e 


(2.20) 


where  Eth,e  denotes  the  electron-initiated  threshold 
energy.  Note  that  ki  is  unknown  and  needs  to  be 
eliminated. 

Using  energy  conservation,  Ei  can  be  expressed  in 
terms  of  the  final  particle  energies  in  the  following  way: 


Ei 


2 2 2 2 2 2 
ft  ki  ft  k2  ft  k3 

“o  + - T.  + - ^ + Err 

2 me  2me  2mh 


2mf 


(2k; 


+ k2  2&) 
x3  mh 


+ Eg 
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2 2 2 
HJSi{2  + iS3_  me  + 
2 me  k?  mh 


(Note : 


JS3 

ki 


mh, 

me 


*2  2 

ft  kj  mh 

2me  me 


+ Eg  — Eth,e 


(2.21) 


Eliminating  ki  from  Eqs . (2.20)  and  (2.21)  leads  to  the 

following  result  for  Eth,e: 

Eth,e  = Eg(l  (2.22) 


Similarly,  one  finds  for  the  hole-ionization  threshold 


Eth,h  — Eg(l 


+ 


mh 

me  + mh 


(2.23) 


From  these  results  some  interesting  conclusions  can  be 
drawn.  As  expected,  the  threshold  energy  depends  primarily 
on  the  bandgap  and  the  effective  masses.  Since  the 
effective  mass  is  a function  of  the  crystal  orientation, 
one  might  expect  that  the  threshold  energy  will  also 
depend  on  crystal  orientation.  Shichijo  and  Hess  [11] 
showed,  however,  that  this  effect  is  lost  in  the 
acceleration  history  of  the  initiating  particle. 

In  addition,  it  is  worth  noting  that  when  me  = mh. 


Eth,e  - Eth,h  - ~ Eg 


(2.24) 


which  is  the  well  known  3/2-  bandgap  rule  [8] . 
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(b)  Three  Simple  Parabolic  Bands 

We  consider  a more  sophisticated  model  by  adding 
another  valence  band  for  light  holes  (shown  in  Fig.  2.2). 
We  assume  that  an  initiating  hole  (ki)  is  located  in  the 
light-hole  valence  band.  After  impact  ionization  we  have 
two  holes  (k2,  k3)  in  the  heavy-hole  valence  band  and  one 
electron  in  the  conduction  band  (ki)  . 

In  terms  of  electron  parameters  we  have  that  before 
the  impact  ionization  process  k2  and  k3  are  occupied. 
After  the  ionization  process  ki  and  ki  become  occupied. 
Hence 


where  mih  and  mhh  denote  the  light-hole  effective  mass  and 
the  heavy-hole  effective  mass,  respectively. 

Before  ionization  the  kinetic  energy  of  the 
initiating  hole  is 


ki  + ki  = k2  + k3 


(2.25) 


which  can  be  rewritten  as 


(2.26) 


2,  2 


2,  2 


(2.27) 


where  Eth,h  represents  the  hole-initiated  threshold 
energy.  After  impact  ionization  we  have  two  holes  and  an 
electron.  From  energy  conservation  Ei  + Ei  = E2  + E3  and 
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(a) 


Figure  2.2  E-k  diagram  with  impact  ionization 
processes : 

(a)  primary  particle  before  ionization 

(b)  resultant  particles  immediately  after 
ionization . 
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with  choosing  the  top  of  the  heavy  hole  band  as  energy 
reference  we  find  for  the  minimum  of  Ei 


Ei 


2 2 2 2 

h ki  ft  k2 

2me  2me 


2 2 
ft  k3 
2mh 


+ E( 


2mhh 


(2k?  + 


2 mhh 
me 


) + Ec 


2 2 2 2 

ft  kvl  mIh  ^vh  _ mlh  _^l_> 

2mih  ^ mhh  k2-  me  v2 


+ Eg 


*2  2 
ft  kyj 

2mih 


<2Yi  i + Y2  % 

*1  *2 


+ Eg 


2 2 

ftkvi  2 
2mih  y1 


) + Eg  - Eth,h 

<2 


(2.29) 


Eliminating  k^  from  Eqs.(2.27)  and  (2.28)  leads  to 


Eth,h 


Eg  ( 1 +~j ) 

y 2 mhh  + me  - mih' 


(2.30) 


Since  the  electron-initiating  ionization  process  is 
identical  to  the  case  of  two  parabolic  bands,  we  have 

Eth,e  = Eg(l  + me  ^mhh)  (2.30) 

(c)  Three  Simple  Parabolic  Bands  with  Spin-orbit 
Splitting  of  the  Light-hole  band 


We  can  make  the  model  more  realistic  by  considering 
two  valence  bands,  one  for  heavy  holes  with  effective  mass 
mh  and  the  other  separated  toward  lower  energy  by  an 
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amount  A due  to  the  spin-orbit  coupling  with  effective 
mass  mih  (Fig.  2.3). 

From  the  momentum  conservation  principle  we  obtain 
the  wave  vector  relationship  where  we  mix  electron  and 
hole  parameters  as  is  usually  done  in  the  literature 


Each  wave  vector  can  be  specified  as  follows: 

ki  : wave  vector  of  initiating  hole  in  the  light  valence 

band 

k2  = k3  = kv,h  (in  the  heavy-hole  valence  band) 
ki  = kc  (in  the  conduction  band) . 

In  addition  one  can  derive  from  Eg.  (2.17) 


Substitution  of  Eqs.  (2.31)  and  (2.32)  into  (2.25)  gives 
ki  = kv,i(“  + — ) (2.33) 


Before  impact  ionization  occurs,  the  total  energy  of  the 
initiating  hole  in  the  valence  band  can  be  expressed  as 


ki  + ki  = k2  + k3 


(2.25) 


(2.32) 


(2.31) 
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Figure  2.3  E-k  diagram  showing  the  impact  ionization 
initiated  by  the  hole  in  the  spin-orbit 
band. 
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2 

Solving  for  k from  Eq.  (2.34)  gives 
j ^2  _ Eth  ~ ^ 2mih 

Vfl  (2/Yi  + 1/Y2)2  t!2 


(2.35) 


From  the  energy  conservation  principle,  we  find  for  the 
minimum  of  Ei 


= Eg  + (K*E) 2-h  + (K-E)i-e 


2 2 2 2 
, iLJSvl  ,9  mih  kvh  _ mh  kc  . 
9 2mih  U mhh  k2  me  k2  } 


ii2kvi , 2 1, 

E9  + 2mih  (7l  + Y2  “ Eth,h 


(2.36) 


In  the  third  line  we  used  Eqs . (2.31)  and  (2.32).  When  we 
substitute  Eq.  (2.35)  into  (2.36),  we  obtain  for  the  spin 
orbit  (s.o),  hole-initiated  threshold  energy: 


Eth,h 


Eg(l 


mih(l  ~ A/Ect) 

2 mhh  + me  - mhh 


(2.37) 


The  electron-initiated  threshold  energy  is  unchanged  and 
has  the  same  form  as  in  the  previous  case: 


Eth,  h 


Eg  (1  + 


me 

me  + mhh 


(2.38) 


For  a material  in  which  Eg  = A,  Eth,h  equals  the 
bandgap,  and  ki  = ki  = k2  = k3  = 0 . In  other  words,  the 
zero  momentum  transferred  in  this  vertical  transition 
reduces  the  threshold  energy  and  increase  p.  This  effect 
is  referred  to  as  resonant  enhancement  of  the  hole- 
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ionization  rate  and  was  experimentally  observed  in 
Al^Gai-xSb  alloys  by  Hildebrand  et  al . [16,17]  near  x = 
0.065.  His  experimental  values  for  p/a  are  plotted  in  Fig. 
2.4.  As  A/Eg  approaches  unity,  (3  becomes  much  larger  than 
a.  For  A/Eg  = 1.02,  p/a  = 20  at  £ = 3 . 3 x 104  V/crrr1.  This 
enhancement  of  p/a  leads  to  a very  low  excess  noise  factor 
for  long  wavelength  (1.55  |im)  APD's. 

Recently,  Brennan  et  al . [18]  found  from  Monte  Carlo 
simulations  that  the  minimum  in  the  hole-ionization  energy 
for  A = Eg  accounts  only  for  a small  part  of  the  resonant 
enhancement  of  p.  Another  important  requirement  for  the 
enhancement  of  P to  occur  is  that  the  hole-ionization 
threshold  energy  (Eth,h)  must  be  very  high  (1.4  eV)  in 
both  the  light-  and  heavy-hole  bands. 

Brennan  et  al . [18]  explained  that  p/a  enhancement 
can  only  be  accomplished  if  enough  holes  are  scattered 
into  the  split-off  band  before  they  ionize  from  either  the 
heavy-  or  the  light-hole  band.  For  that  to  occur,  it  is 
necessary  for  the  holes  to  drift  to  a sufficiently  high 
energy  where  the  density  of  states  of  the  split-off  band 
is  appreciable.  Therefore,  in  order  for  the  resonance  to 
occur,  Eth,h  must  be  high  for  both  the  heavy-  and  light- 
hole  bands  (~  1.4  eV) . 


0 / a 
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Figure  2.4  The  ionization  coefficient  ratio  p/a 
plotted  versus  A/Eg  (Ref.  17) 
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(d)  A Nonparabolic  Conduction  Band  and  One-parabolic 
Valence  Band 

The  conduction  band  central  valley  in  III-V  materials 
can  be  quite  well  represented  by 


2 2 
-h  k 

2me 


= E (1  + a'E) 


(2.39) 


where  a'  is  a measure  for  the  deviation  of  the  E(k)  band 
from  a parabolic  shape.  Ridley  [19]  showed  that  the 
threshold  energy  for  electrons  becomes  in  this  case 


Eth,  e 


Eg( 


1 + n ; 


[1  + 


M-d  + 2H) 
(1  + ^i)2 


a'Eg] 


(2.40) 


where  |i  is  me/mh,  and  me  and  mh  are  the  electron  and  hole 
effective  masses  respectively. 

For  example,  for  GaAs , fi  = 0.14  9 and  a = 0.64  eV-1. 
Substitution  of  these  values  into  Eq.  (2.40)  yields  Eth,e 
= 1.82  eV . Landsberg  and  Yu  [20]  found  a more  accurate 
expression  for  the  threshold  energy  of  electron-electron- 
hole  impact  ionization,  given  by 


Eth,  e ~ Eg- 


UL 

2\\a' 


4ll(l+2n) 

(l+)l)2 


a'Eg  ] 


(2.41) 


When  we  substitute  above  values  into  Eq.  (2.41),  we  obtain 
Eth,e  = 1.95  eV,  which  is  in  better  agreement  with 
experimental  findings.  The  derivation  and  a more  detailed 
discussion  of  expression  (2.41)  will  be  presented  in 
Chapter  IV. 
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2.2  Noise  Theory 

2.2.1.  The  Binomial  Process 

Let  an  experiment  be  tried  m times.  Let  p be  the 
probability  that  the  trial  is  successful  and  1 - p the 
probability  that  it  fails.  Let  there  be  n successes  for 
the  m trials,  then 

n = mp  (2.42) 

var  n = mp (1  - p) 

and  the  probability  distribution  is  given  by 
_ , . m! 

Pm(n)  = n7  (m  - n)  ! pn(1  ~ p>m_n  (2.43) 

This  is  called  the  binomial  law. 

For  example,  in  terms  of  a photodiode,  a "trial"  is 
equal  to  the  absorption  of  a photon  by  the  p-n  junction, 
"success"  is  when  the  generated  electron-hole  pair  is 
collected. 

2.2.2.  Poisson  Process 

The  photon  detection  by  a photodiode  or  the  emission 
of  a photon  in  the  visible  wavelength  spectrum  can  be 
described  as  a Poisson  process  [21] . If  the  arrival  of 
photons  at  the  detector  can  be  described  as  a series  of 
independent  events  occurring  at  random  instants,  then  the 
binomial  distribution  reduces  to  the  Poisson  distribution 
as  will  be  shown  in  the  following. 
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If  m goes  to  infinity  and  p goes  to  zero  (mp  = n is 
kept  constant),  then  Eg.  (2.43)  can  be  approximated  as 


Pm(n)= 


m (m-1) • • • (m-n+1) 


n! 


mn  pn  ( 1 -p ) m~n 
n! 


n 


a-p)m'n 


pn(l-p)m  n (where  we  divided 
out  the  term  (m-n) ! ) 


Define  (1-p) m-n  = X/  then  (m-n) In ( 1-p) =ln  (x) , or  with  p -» 
0 In  (1-p)  ~ -p  and  x = e-mPe+nP  = e-n-l.  Consequently,  we 
derive  at 


P(n) 


in ) 
n! 


n 

-^exp  (-n) 


(2.44) 


m! 

n!  (m  - n)  ! 


Pn(l 


- p)  m-n  . e-mBi2Eli 


(2.44) 


which  is  the  Poisson  distribution.  The  variance  of  a 
Poisson  distribution  is  given  by 

var  n = n (2.4  6) 

2.2.3.  Variance  Theorem 


We  consider  the  electron  multiplication  process  as  it 
occurs  in  the  avalanche  region  of  an  APD. 

Let  N be  the  number  of  electrons  arriving  at  the 
device.  Let  Xi  be  the  number  of  secondary  electrons 
(offspring)  produced  the  ith  primary  electron.  Let  a 
second  quantity  n be  defined  by 
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N 

n = 2 xi  (2.47) 

i=l 

Since  all  of  the  primary  electrons  have  the  same  a priori 
probability  for  multiplication,  it  follows  that 

Xi  = M ; X?  = x2  (2.48) 

independent  of  i,  where  M is  the  gain.  Let  N and  Xi  both 
fluctuate;  then 

n = NX  = NM  (2.49) 

var  n = M2  var  N + N var  X 

which  is  called  Burgess'  variance  theorem  [21]. 

If  we  make  use  of  the  fact  that  var  N = "n  (Poisson), 
then  Eq.  (2.48)  reduces  to 

var  n = NX2  (2.50) 

Assume  that  the  currents  flowing  in  the  device  show  full- 
shot  noise,  then  we  have  the  following  correspondences: 

var  n corresponds  to  Sid  = 2qld:  Sid  is  the  noise 

spectral  density 
of  diode  current 

var  N corresponds  to  Sipr  = 2qIpr:Sipr  is  the  noise 

spectral  density 
of  primary  current . 

Then  Eq.  (2.48)  can  be  rewritten  as 

Sld  = m2  SIpr  + SIpr  var  x (2.51a) 

or 

Sid  = sIpr  <m2+  var  x)  (2.51b) 


The  first  term  of  Eq.  (2.51a)  represents  amplified  primary 
noise,  and  the  second  term  is  identified  as  multiplication 
noise  [21] . 

Tager  [22]  and  McIntyre  [23,24]  showed  that  var  X is 
given  by 

var  X = M2  (M-l)  : Tager  (equal  ionization  probabilities) 

(2.52) 

var  X = M(M-l)  + k(M-l)2M:  McIntyre  (2.53) 

where  k = p/a 

If  we  substitute  Eq.  (2.52)  into  (2-51),  we  get 

SId  = sIpr  [m2+  M2  (M-l)]  = SiprM2  (2.54) 

This  is  the  result  Tager  obtained  [22]  . In  the  following 
sections  Eqs.  (2.52)  and  (2.53)  will  be  discussed  in  more 
detail . 

2.2.4.  Excess  Noise  Factor.  F 

We  rewrite  Eq.  (2.51)  by  using  the  definition  of  a 


variance  and  obtain 
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where  F is  known  as  the  excess  noise  factor.  F is  a 
measure  for  the  amount  by  which  the  noise  exceeds  that  of 
an  ideal  current  multiplier.  To  clarify  the  meaning  of  F, 
we  will  consider  the  following  example. 

Let  an  input  signal  have  the  current  value  Iph;  then, 
after  multiplication  the  output  signal  current  becomes 

lout  = IphM  (2.56) 

If  we  have  an  ideal,  noise-free  amplifier,  the  input 
current  noise  would  be  amplified  by  a factor  M;  hence 

V 2 qi^Af  output  — V^qlphAf-M  (2.57) 

Squaring  both  terms  in  Eg.  (2.57)  yields 


slph  = siprM2  (2.58) 

By  comparing  Eq.  (2.58)  with  (2.55),  we  can  deduce  that 
for  an  ideal  current  multiplier  the  excess  noise  factor 
reduces  to  unity. 

However,  in  reality  nonideality  exists  for  which  the 
second  term  (multiplication  noise)  of  Eq.  (2.50)  is 
responsible.  Therefore,  the  noise  level  increases  by 
(MF) 1 /2 . 

The  excess  noise  factor  F is  always  greater  than 
unity  and  is  an  increasing  function  of  M.  The  value  of  F 
depends  upon  the  detailed  characteristics  of  the  avalanche 
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junction  and  upon  the  nature  of  the  primary  excitation: 
hole  injection,  electron  injection  or  mixed  injection. 

Tager  [22]  showed  that  for  equal  ionization 
coefficients  of  electrons  and  holes  the  noise  spectral 
density  is  given  by 

Sid  = 2qIprM3  (2.53) 

where  it  is  assumed  that  the  incoming  particles  show  full- 
shot  noise.  By  equating  Eqs . (2.53)  and  (2.55)  we  have  F 

=M. 

McIntyre  extended  the  analysis  to  the  case  of  unequal 
ionization  coefficients  [23] . He  assumed  that  the  relation 
between  these  two  is  k = (3/a  and  is  independent  of  the 
field  strength.  The  noise  spectral  density  is  then  given 
by 

Sid  = 2qIprM3[l  - (1-k)  (^)2]  (2.60) 

By  equating  Eqs.  (2.60)  and  (2.55)  we  can  obtain  the 
expression  for  the  excess  noise  factor 

(M-l) 2 

F = M [ 1 - (1-k)  m2;  ] (2.61a) 

whereas  Tager' s theory  results  in 

F = M (2.61b) 

For  later  use  we  rewrite  the  excess  noise  factor  as 


follows : 
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(2.62) 


(2.63) 


Eq.  (2.62)  follows  from  Eq.  (2.55),  and  Eq.  (2.63)  is 
obtained  by  equating  Eq.  (2.51)  and  Eq.  (2.55). 

2.2.5.  Generating  Functions  and  Compounding  Theorem 
1)  Generating  Functions 

For  random  variables  whose  only  possible  values  are 
non-negative  integers,  for  example,  0 and  1 in  the 
Bernouilli  process  and  0,  1,  2,  • • • in  the  binomial 

counting  process,  a function  related  to  the  characteristic 
function  is  the  generating  function  [25]  of  the  Fowler- 
type  of  variable  k defined  by 


here  pk  = probability. 

The  generating  function  of  a non-negative,  integer 
valued,  random  variable  k is  formally  related  to  the 
characteristic  function  <J>  of  k through  the  change  of 
variable  s = eit:. 


oo 


g(s)  = X Pksk 
k=0 

= E[sk] 


(2.90) 


<J>(t)  = E[ei1:] 


= E[(eit)x] 

= g(ei1:) 


(2.91) 
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Thus  generating  functions  have  the  three  basic  properties 
of  characteristic  functions: 

a)  A generating  function  determines  the 
distribution  function  uniquely. 

b)  The  generating  function  of  a sum  of 
independent,  non-negative,  integer  valued,  random 
variables  is  the  product  of  their  generating 
functions . 

c)  The  moments  may  be  obtained  through 
successive  differentiation. 

2)  Compounding  Theorem 

Let  N,  Xi,  X2,  • • • be  independent  non-negative 

integer  valued  random  variables  and  suppose  we  wish  to 
determine  the  generating  function  gR(s)  of  the  sum  R = Xi 
+ X2  + • • • + Xn  . 

Let  gu(s)  be  the  generating  function  of  N and  suppose 
that  all  Xi  have  the  same  distribution  functions  with  the 
common  generating  function  g(s). 

Then, 

gR(s)  = E [sR] 

= E[sx1+x2+'  • ,+xN] 

= E{E[sx1+x2+‘  ■ •+XN|N|} 

oo 

= Xe[sx1+x2+-  • ■ +XN  | N=n ] Pr  { N=n } 
n=0 
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oo 


(since  N and  Xi  are  independent) 


oo 


= X g(s)nPr{N=n} 
n=0 


= E [g  (s) N] 


(2.92) 


Therefore, 


gR(s)  = gN[g(s>] 


(2.93) 


This  is  the  compounding  theorem.  If  we  consider  the 
avalanche  multiplication  process,  then 

X corresponds  to  the  number  of  offsprings, 

N corresponds  to  the  number  of  incoming  carriers. 


R corresponds  to  the  total  number  of  carriers,  M, 
collected  after  multiplication. 

The  compounding  theorem  is  a very  useful  tool  to 
analyze  the  branching  process  in  APD ' s and  is  used  as  a 
central  tool  by  C.M.  Van  Vliet  et  al.  [26,27]  to  formulate 
the  multiplication  noise. 

2.2.6.  General  Theory  of  Avalanche  Noise  by  Van  Vliet  and 


Contrary  to  the  standard  theories  [22,23,24,28]  which 
treat  the  impact  ionization  of  carriers  as  a continuous 
spatial  process.  Van  Vliet  et  al.  [26,27]  developed  a new 
general  theory  employing  the  method  of  recurrent 
generating  functions.  This  theory  is  quite  general,  since 


and 


Rucker 
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it  covers  the  single-carrier  process  as  well  as  the  two- 
carrier  process,  from  N = 1 to  N = °°,  where  N is  the 
number  of  possible  multiplications.  Therefore,  the  results 
by  Tager  [22],  McIntyre  [26,27],  Personick  [28],  and 
others  are  found  to  be  special  cases  of  this  theory.  One 
of  the  important  features  of  Van  Vliet's  theory  is  that  it 
properly  describes  the  transition  from  continuum  device 
physics  to  discrete  device  physics.  The  description  of 
this  transition  is  important  in  many  modern  devices.  Since 
this  theory  is  based  on  a finite  number  of  discrete 
ionization  processes,  the  result  is  also  entirely 
discrete.  Van  Vliet  predicted  the  presence  of  discrete 
steps  in  the  noise  for  a finite  number  of  possible 
multiplications  in  the  avalanche  region.  J.  Gong  [29] 
verified  this  by  measurements. 

In  this  section  we  will  outline  the  principles  and 
the  results.  The  theory  starts  from  the  following 
assumptions : 

1)  Ionization  can  take  place  after  an  energy  gaining 
path  li,  with  li  = Eth/ql^il/  where  Eth  is  the  ionization 
threshold  energy  and  Ei  is  the  average  electric  field  in 
the  interval  li. 

2)  The  number  of  possible  ionizations  per  carrier 

transit  is  finite  since  we  can  divide  the  avalanche  region 

N 

w into  N sections  of  length  li,  such  that  w = £li. 

i=l 
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3)  The  a priori  chance  for  electron-initiated 
ionization  after  an  energy-gaining  path  li  is  X-i.  In  the 
same  fashion  we  define  a priori  chance  fig  for  holes. 

4)  After  each  collision  all  energy  is  lost  and  the 
momentum  is  randomized,  so  that  each  possible  ionization 
during  the  successive  paths  li  • • • In  is  governed  by  an 
independent  Bernouilli  trial.  Note:  This  assumption  needs 
to  be  modified  a little  bit  according  to  the  theories 
[11,12]  discussed  in  the  previous  section  (see  impact 
ionization  mechanism) . 

A.  One-carrier  Process  Let  N be  the  number  of 
possible  ionizations  per  primary  carrier  transit  and 
consider  only  the  electrons  in  this  process  as  shown  in 
Fig.  2.5.  We  consider  the  subtrees  arising  from  the  first, 
second,  • • •,  Nth  ionization  event  of  the  incoming 
carrier.  These  subtrees  are  denoted  by  Yn,  Yn-i,  • • •,  Yi, 
respectively;  thus 

X = XN  = 1 + Yn  + Yn-i  + • • • + Yi. 

We  now  observe  that  the  subtree  Yn  is  identical  to  the 
tree  X,  except  that  only  a maximum  of  N - 1 ionizations 

can  occur  and  that  the  subtree  is  "fed"  by  a Bernouilli 
trial  (Fig.  2.6)  . If  gN  and  <))pg  denote  the  generating 

functions  of  Yn  and  Xn,  respectively,  then  the  compounding 
theorem  gives  [see  Appendix  A] 


gN  (z)  = co [ <J>N — i (z)  ] 


(2.94) 
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Figure  2.5  Impact  ionization  process  involving 

N possible  ionizations  per  transit  of 
the  primary  carrier. 
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Sub-Tree 


Sub-Tree 


Sub-Tree 


Sub-Tree 


XN  = 1 + XYi 


i = 1 


Figure  2.6  Diagram  of  the  one-carrier  multiplication 
process . 
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where  co  is  the  generating  function  for  a Bernouilli  trial 
with  a priori  chance  X± . Thus 


In  general,  for  branch  i 

gN— i+i(z)  = C0i[(|)N_i(z)]  = X,i<|)N_i  ( z ) - Ad  + 1(2.96) 

We  also  know  from  looking  at  Fig.  2.6, 

N 

XN  = 1 + £ Yj 

j=l 

Then  in  terms  of  generating  functions 


gN(z)  - M>N-1  “ A,  + 1 


(2.95) 


N 


<t>N  = z X 

j=l 


(2.97) 


where  gQ  = z . In  a similar  way 


N— l 


^n-i  = 2 X gj 
j=i 


(2.98) 


Dividing  Eg.  (2.97)  by  (2.98),  we  have 


<1>N 


(2 . 99) 


By  using  Eqs . (2.95)  and  (2.99)  we  obtain 


<I>N  (Z ) = (Z)  [Axt>N_i  (z)  - X + 1]  (2.100) 


which  is  solved  with  the  boundary  condition  <])o(z)  = z.  The 
following  result  for  the  multiplication  factor  M is  found 
by  iterating  Eq.  (2.100), 
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M = <^(1)  = (1  + X)N  (2.101) 

For  the  noise  we  find 

var  X = (t^(l)  -[<to(l)]2  + <^(l)  = M(M-1} 

(2.102) 

Substitution  of  Eq.  (2.99)  yields 

var  XN  = (2  - m1/N)m1~1/N(m  - 1)  (2.103) 

For  N = 1,  this  becomes 

var  Xi  = (2  - N) (M  - 1)  (2.104) 

For  N — >00, 

var  Xoo  = M(M  - 1)  (2.105) 

which  is  the  Tager's  result.  The  multiplication  factor  can 
also  be  found  for  the  asymptotic  case  N — > 00.  since  X = 
(Xw/N,  where  w/N  = 1,  Eq.  (2.99)  reduces  to 

Moo  = lim  (1  + ^)N  = eaw  (2.106) 

N— W 

TwQ~carrigr Process . First  we  assume  that  a 

primary  electron  comes  in  at  a position  (k,l)  where  it  can 
create  a hole-electron  pair.  Here  k is  the  number  of 
possible  subsequent  ionizations  by  the  electron  which 
moves  to  the  right  and  1 is  the  number  of  possible 
ionizations  by  the  hole  which  moves  to  the  left,  thus  k + 
1 — N,  see  Fig.  2.7.  The  total  number  of  hole-electron 
pairs  due  to  the  pair  coming  in  at  (k,l),  including  this 
original  pair,  will  be  denoted  as  Xkfi  and  given  by 


hoton 
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Figure  2.7  Diagram  of  the  two-carrier  multiplication 
process . 
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Xk,l  - 1 + Yk-1,1+1  + Yk-2,1+2  + * ’ • + y0fN 

+ Vk+1,1-1  + Vk+2,1-2  + • • • + yNf  o (2.107) 

Let  ^k, l ( z ) he  the  generating  function  of  Xk,  then  (|)k  x(z) 
can  be  expressed  by  the  product  of  generating  functions  of 
subtrees  as  follows: 

^k,  1 = zgk-l , 1+1 ' gk-2 , 1+2 " ‘ ' ‘9k+l,  1-1 '9k+2,l-2'  ’ '9n,0 

(2.108) 

or 


k 1 


^i^  z n^-j.i  j n^j.i-j 

(2.109) 

j=i  j=i 

By  letting  k — > k + 1,  1— >1-1, 

k 1 

^k+1,1-1  (z)  =z  J'J^k+l-j  , l-l+j  ']"~J<?k++l+j  , l-i-j  (2  . HO) 
j=l  j=l 


Dividing  Eg.  (2.109)  by  (2.110),  we  find  the  ratio 


^k , 1 gk+l,l-l 

^k-1,1+1  gk,l 


(2.111) 


Now  the  generating  function  gk_1/1+1  of  yk-i,i+i  is 
identical  to  that  of  Xk_-|^  1+1  except  that  the  "feeding"  is 

by  a Bernouilli  trial  of  a priori  chance  Xi+i;  thus 
according  to  the  compounding  theorem,  it  is  given  by 
^l+l^k-i,  1+1  - A.i+1  + 1.  likewise,  the  subtree  yk+i  i_1  is 


45 


identical  to  that  of  Xk+1^1_1,  except  that  the  feeding  is 

by  a Bernouilli  trial  of  a priori  chance  |li-i,  resulting 
in  a generating  function  i_]_  “M-i-i  + i.  Therefore, 

Eq.  (2.111)  reduces  to  the  recurrent  relation  as 

^k,  1 l(z)  ^1  + 1]  - ^+1, 1-1  [^l-l^k+l,  1-1  " M-l-1  + 1] 

(2.112) 

When  k = N,  1=0,  we  have  the  case  of  an  electron-hole 
pair  coming  in  at  the  left.  However,  the  hole  so 
introduced  is  immediately  carried  off  at  the  negative 
electrode  at  the  left.  So  this  situation  is  equivalent  to 
the  case  that  only  an  electron  comes  in  at  the  left  and 
collected  at  the  right  after  multiplications.  Therefore, 
once  we  know  the  generating  function  <])N,0r  we  can  find  the 
mean  (i.e.,  multiplication  factor,  M)  and  var  X by 
calculating  the  first  and  second  derivatives  of  <|)n,o(z)  at 
z = 1,  respectively.  Thus 

M = <Xn,0>  = , o ( 1 ) (2.113) 

var  XN,0  - 4>n,  0 (!)  ~ + ^(D 

= <1>n,0  (D  - M2  + M (2.114) 

Eq.  (2.112)  is  used  to  find  the  generating  function  <])n,0- 
But  the  algebra  involved  in  Eqs . (2.113)  and  (2.114)  is 

quite  long  and  complex.  Hence  we  will  skip  the  derivation 
and  show  the  important  results  below. 
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M = 


(1+X) N (l-k) 


(l+kX)N+i  - k(l+A.)N+i 


where  k = ^ (2 . 115) 


Moo  =■ 


a - 3 


a - 3 


1 + 3w/N  M _ o g exp(3-g)  - 3 
^1  + gw/N  p 


(2.116) 


..  -t.1  « 1— kX^  1+X  1 

^arX-" {-X  + 2 -TTi-[Mnk  — + } (2.117) 


Efr  (Mn-D  (l-k) 

2 + X + kX 


1+kX 


From  Eq.  (2.117)  we  can  derive  the  following  limiting 
cases : 

/ s ^ t X ( 1— X+3|i+Xji.) 

(a)  N - 1,  var  X = — — 7t7o~  (Lukaszek)  (2.118) 


(1-Afi) 


(b)  N = °°,  var  Xoo  = M(M-l)  + kM(M-l) 2 (McIntyre) 

(2.119) 

1 _ x 

(c)  |i  = 0,  var  X = y~ ~ ^ M(M-l)  (one-carrier  process) 

(2.120) 

2.2.1.  Discrete  Noise  Characteristics  of  the  Avalanche 
Photodiode 

Many  theories  of  avalanche  noise  (Tager  [22]  McIntyre 
[23]  and  Personick  [28])  only  deal  with  situations  where 
the  number  of  possible  ionizations  is  very  large.  However, 
in  many  modern  devices,  where  the  depletion  region  is 
thin,  the  number  of  possible  ionizations  is  finite  and 
very  small  (N  = 1 ~ 5)  . 

As  shown  in  the  previous  section.  Van  Vliet  and 
Rucker  developed  a general  theory  which  can  be  applied  to 
such  cases.  To  visualize  the  difference  between  two  cases. 
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we  need  to  plot  the  curve  for  var  X versus  multiplication 
factor  M. 

Eq.  (2.117)  gives  the  expression  for  var  X.  To  find 
var  X as  a function  of  Mn,  we  must  eliminate  X with  the 
help  of  Eq.  (2.115).  The  results  so  obtained  are  shown  in 
Fig.  2.8.  It  is  clear  that  the  noise  is  at  all  times  below 
the  McIntyre  limit.  It  is  also  possible  to  construct  an 
overall  curve  for  var  X in  which  N is  the  dependent 
variable. 

The  idea  is  as  follows.  For  a given  N,  say  N = 1,  when 
the  field  increases  A,(l)  increases,  and  so  does  var  Xn- 
Since  X is  bounded,  0 < X < A.max  < 1,  var  X cannot 
increase  without  switching  to  N = 2 if  X arrives  at  ^max. 
Thus,  the  process  switches  over  to  N = 2,  with  X{2)  being 
less  than  A,i,max  in  order  to  maintain  the  same  gain;  then 
X{2)  increases  up  to  A-2,max,  etc.  Therefore,  the  plot  of 
var  X vs.  M shows  break  points  at  such  M values  where  X is 
discontinuous,  due  to  the  regime  switching  from  N to  N + 
1.  Fig.  2.8  is  obtained  by  piecewise  fitting  the  results 
for  N = 1 , 2,  ■ • • etc. 

In  1980  J.  Gong  [29]  measured  the  APD  noise 
characteristics  and  experimentally  verified  the  above 
theory.  Instead  of  var  X,  he  made  a plot  of  var  X + vs. 
M,  which  is  identical  to  a plot  of  ieq  vs.  Id.  To  see 
this,  we  recall  Eq.  (2.63) 
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Figure  2.8  A plot  of  the  var  X versus  multiplication 
factor  M (Ref.  29) . 
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F 


1 + 


var  X 
M2  ' 


This  equation  can  be  rewritten  as 


FM2  = var  X + M2  = x2 

If  the  primary  current  has  shot  noise,  then 

Sip  = 2qIprX  = 2qleq 

From  this  equation,  we  obtain  the  ratio 


(2.121) 


Also 


t22,  = var  X + M2 . 
-Lpr 

Id 


Lpr 


M 


(2.122) 


(2.123) 


Dividing  Eq.  (2.112)  by  (2.113),  we  have 


lea  _ var  X + M2 
Id  M 


(2.124) 


Therefore,  plotting  var  X + M2  vs.  M is  a more  convenient 
way  to  characterize  the  APD  noise.  One  of  Gong's  plots  is 
shown  in  Fig.  2.9  as  an  example. 

2.2.8.  Excess  Noise  Factors  for  the  Staircase  Avalanche 
Photodiode 

The  staircase  APD  (SAPD)  was  proposed  by  Capasso  [30] 
and  considered  to  be  a solid-state  analogy  of 
photomultiplier  tubes  (PMT's).  A detailed  description  of 
this  device  will  be  given  in  section  III.,  1-3.  In  this 
section,  we  are  interested  only  in  the  derivation  of  the 
excess  noise  factor  F of  the  SAPD. 
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Figure  2.9  A plot  of  I 


eg  versus  Id  (Ref.  29) . 
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a)  Derivation  of  F by  Capasso.  The  band  structure  of 
a SAPD  is  shown  in  Fig.  2.10.  Consider  a photoelectron 
generated  next  to  the  p+  contact.  Under  the  combined 
action  of  the  bias  field  £ and  the  grading  field  AEc/ql, 
where  Aec  is  the  conduction  band  step  and  1 is  the  width 

of  one  step,  the  electron  drifts  toward  the  first 

A£  o 

conduction  band  step.  The  combined  field  £ - —p  is  small 

enough  so  that  the  electron  does  not  impact  ionize  before 
it  reaches  the  step.  After  the  step,  since  Aec  = Eth,e^ 
the  electron  creates  a hole-electron  pair.  This  ballistic 
ionization  process  is  repeated  at  each  step. 

Ideally,  the  avalanche  gain  per  stage  is  2.  In 
practice,  the  gain  is  2 - 9,  where  d is  the  fraction  of 
electrons  which  do  not  impact  ionize.  If  N is  the  number 
of  stages,  the  total  gain  is  given  by 

<M>  = (2  - 9)N.  (2.125) 

When  we  neglect  the  dark  current,  the  noise  spectral 
density  is 

Sid  = 2qIpr<M>2F(N,a)  , (2.126) 

where  F is  the  excess  noise  factor  and  is  a function  of  N 
and  d.  If  we  write  the  shot  noise  generated  by  the  nt*1 
step  as 


Si, n = 2qIn_!var  m, 


(2.127) 
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Figure  2.10  Band  digram  of  the  staircase  avalanche 
photodiode  (SAPD) (Ref.  6) : 

(a)  SAPD  is  unbiased. 

(b)  SAPD  is  reverse  biased. 
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where  m is  the  gain  of  the  nth  stage.  Since  we  will  deal 
with  a binomial  process  ( see  section  2.2.1.), 

<a2>  = a 

<m2>  = <(2-3)2>  = <4  - 43  + 32>  = 4-33  (2.128) 

<m>2  = (2-3)2  = 4 - 43  + 32  (2.129) 

Using  Eqs.  (2.128)  and  (2.129),  eq  (2.127)  becomes 

Si,n  = 0 “ 52)  (2qln-i)  (2.130) 

where  In-i  = lo(2-3)n_1  is  the  average  current  into  the  nth 
stage . 

We  assume  that  the  noise  contribution  of  the  n^h 
stage  (1  < n < N)  is  multiplied  noiselessly  by  the  rest  of 
the  stages  by  a factor  M2,  where  M2  is  given  by 

M2  = [ (2-3)  N-nj  2 = (2-3)  2 (N-n)  (2.131) 

The  total  output  noise  per  unit  bandwidth  is  then 

Sid  = 2ql0(2  - 3)2N  + Si,  i(2  - 3)2(N-1) 

+ Si,  2 (2  - 3)2  (N-2) 

+ . . . 

+ Si, n (2  - 3)  0 

N 

= 2ql0  (2  - 3)2N  + X2qIo  (2-3)  n-l(3-32)  (2-3)2  (N-l) 

n=l 
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N 


2ql0  (2  - 8)2Htl  + ^(2_13)d  (2.132) 


n=l 


By  using  the  following  summation  formula  for  this 
geometrical  series 


SxT  “ 


al  ~ raN 


jn  - 1 - r 

Eq.  (2.132)  reduces  to 


r < 1 


(2.133) 


Sid  - 2ql0 (2  - a)2N[i  + d{1  2 ^ N}j.  (2.134) 


Therefore,  the  excess  noise  factor  becomes 


F(N,3)  = 1 + 


d{l  - (2-5) -N} 

2“^ 


(2.135) 


This  result  is  plotted  in  Fig.  2.11.  If  5 = 0,  i.e.,  all 
the  electrons  ionize  at  each  step,  F becomes  unity. 
Therefore,  this  SAPD  can  in  principle  be  a noise-free 
device . 

*>)  • Application  of  the  Theory  bv  Van  Vliet  et  al . to 
SAPD ' $ Van  Vliet  et  al.  [26,  27]  assumed  that  both 
electron-  and  hole-impact  ionization  occurs  at  discrete 
positions  (see  section  II.,  2-7)  . However,  in  SAPD's 
electrons  impact-ionize  at  each  step  while  holes  impact- 
ionize  at  any  position  in  the  depletion  region.  Therefore, 
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MULTIPLICATION  PER  STAGE 


Figure  2.11  Excess-noise  factor  of  the  SAPD  as  a 

function  of  the  multiplication  per  stage 
N represents  the  number  of  staaes 
(Ref.  6) . y ' 
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we  need  to  modify  the  assumption  of  impact  ionization 
rates  as  follows. 

X:  electron-impact  ionization  rate  per  stage 
|i:  hole-impact  ionization  rate  per  stage. 

A,,  (j.  are  given  by 

k = 1 - d (2.136) 

and 

L 

[i  = exp [ Jpdx]  - 1 (2.137a) 

0 

respectively,  where  L is  the  width  of  each  step. 

If  P is  very  small,  then 

M-  = (2.137b) 

where 

L 

F = l (2.138) 

The  derivation  of  Eq.  (2.137a)  goes  as  follows.  We  assume 
that  impact  ionization  caused  by  electrons  between  the 
steps  is  negligible.  Then  we  can  use  the  treatment  of  one- 
carrier  process  associated  with  holes.  The  hole  current 
rate  equation  can  be  expressed  by 

dip  (x)  = plpdx.  (2.139) 

We  rewrite  this  as 


Th  solution  of  this  equation  is  simply  given  by 


(2.140) 
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L 

Ip(x)  = Ip  (0)  exp(Jpdx')  (2.141) 

o 

To  find  the  multiplication  factor  for  one-stage,  we  divide 
the  output  current  at  location  L by  the  input  current,  and 
we  have 

L 

Mp  = Ymf  = exp[  Jpdx'  ] = ePL  (2.142) 

v ; o 

where  T is  given  by  Eq.  (2.138).  Since  an  one-carrier 
process  is  involved,  we  obtain  for  the  expression  of  the 
multiplication  factor  for  one-stage  in  terms  of  the  hole- 
impact  ionization  rate  per  stage  using  Eq.  (2.101) 

Mp  = 1 + \i  (2.143) 

From  Eqs . (2.142)  and  (2.143)  we  obtain  the  result  as 

L 

|l  = exP [ Jpdx]  -1  (2.137a) 

o 

If  we  expand  the  exponential  terms  in  series  form, 

Mp  = ePL  = 1 + “PTL  + (~^)2  + • • • (2.144) 

For  'P’L  « 1,  we  have 

Mp  = 1 + ITl.  (2.145) 

From  Eqs.  (2.143)  and  (2.145)  we  also  obtain  the 
approximate  result 


H = 'PL 


(2.137b) 
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If  we  neglect  hole-impact  ionization,  then  this 
situation  becomes  a one-carrier  process.  Therefore,  Eq. 
(2.101)  gives  the  gain  of  SAPD  as 
M = (1  + A)N 

Using  Eq.  (2.136),  this  reduces  to 

M = (2  - 3)N  (2.146) 

which  is  the  same  result  obtained  in  the  previous  section. 

For  the  excess  noise  factor  F we  find,  using  Eqs . 
(2.101)  and  (2.102), 


F = 1 + 


var  X 


M2 

„ 1 - A,  M 

= 1 + TTT  m2  (M  - 


(2.147) 


Substituting  of  Eq.  (2.136)  and  Eq.  (2.146)  into  (2.147) 
yields 

f = i + 2-5  [i  - (2  - arN] 

which  is  the  same  result  derived  by  Capasso  in  the 
previous  section.  This  derivation  was  published  by  Teich 
et  al . [31] . 

2.3.  DEVICE  OPERATION 
2.3.1.  Design  Considerations 
2. 3. 1.1  Gain  (Multiplication) 

Assume  that  an  electron  current  ie(0)  is  injected 
into  a depletion  layer  of  width  w at  the  plane  x = 0.  We 
will  assume  that  the  field  strength  is  high  enough  to 
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produce  an  avalanche  and  that  there  is  negligible  thermal 
or  photo-generation  of  carriers  in  this  region.  At  any 
point  x such  that  0 < x < w,  the  rate  of  generation  of 
carriers  is  given  by 

die (x) 

—  = aie(x)  + pih(x)  (2.148) 

Also  at  any  point 

ie  (x)  + ih(x)  = I = constant  (2.149) 

Eqs . (2.148)  and  (2.149)  give 
die  (x) 

—  (a-p)ie(x)  = pi  (2.150) 

Eq.  (2.150)  is  of  a standard  form  which  integrates  to  give 

X X 

ie(0)  + J Pi  exp[-J  (a-p)dx'Jdx 

ie(x)  2 (2.151) 

X 

exp[-J  (a-P)dx'] 

0 

Now  we  define  the  multiplication  factor  for  the  injected 
electrons  as 

W X 

ie(0)  + J Pi  exp[-J  (a-p)dx']dx 

,,  _ 1 _ ^e(w)  0 0 

ie  (0)  ie  (0)  w 

exp[-J  (oc-P)dx'] 

0 


By  rearranging,  we  get 


(2.152) 
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Mn 


1 

w x 

1 - J a exp [- 1 (a-p)dx']dx 
0 o 


(2.153) 


In  a similar  way,  we  determine  the  multiplication 
factor  or  gain  for  injected  holes  as 


exp[-J  (a-p)dx] 

MP  w ^ (2.154) 

1 - J a exp  [ - J (a-p)dx']dx 
0 o 

Eqs . (2.153)  and  (2.154)  will  be  used  to  determine  the 

ionization  coefficients  a and  p as  explained  in  the 
following  sections. 

2. 3. 1.2.  Noise  Equivalent  Power  and  Detectivity 

Consider  a detector  receiving  modulated  radiation. 
Let  a radiant  power  p produce  a reading  R of  the 
instrument.  Then  we  define  the  technical  sensitivity  by 


[21], 


C 


R 

P 


If  a detector  produces  a voltage  V,  then 


(2.155) 


where  the  unit  of  C is  volt/watt.  If  a detector  produces 
current  I,  then 


C = CT  =j 


(2.156) 


61 


To  define  the  noise  equivalent  power  (NEP) , we  assume 
that  the  noise  of  the  detector  plus  indicating  instrument 
has  a spectral  intensity  SR(f).  Then,  the  NEP  Peq  is 
defined  by 


Peq 


[ Sr ( f ) ] 1/2 
C 


(2.157) 


where  C is  given  by  Eq.  (2.155)  or  Eq.  (2.156).  The  NEP 
Peq  is  expressed  in  watt/VHz,  and  is  the  power  that  gives 
the  same  reading  as  the  noise  per  unit  bandwidth. 

Usually,  Peq  is  proportional  to  the  square  root  of 
the  detector  area  A.  We  can  normalize  the  value  of  Peq  to 
a unit  area  by  setting 


P*  _ peg 
a1/2 


(2.158) 


where  P@q  is  expressed  in  watt/  (Hz1/2cm)  . 

The  best  detector  type  is  the  one  with  the  smallest 
Peq.  Since  one  wants  to  identify  a good  detector  by  a 
large  number,  we  define  the  detectivity  D*  by  [21] 


D* 


ai/2 


eq 


(2.159) 


which  is  expressed  by  (Hz1/2cm) /watt . 

2. 3. 1.3.  Bandwidth 

The  response  time  of  an  APD  is  determined  primarily 
by  four  effects  [32] : a)  the  avalanche  build-up  time,  b) 
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the  drift  time  of  the  carriers  across  the  depletion  width, 
c)  the  diffusion  time  of  the  generated  minority  carriers 
towards  the  active  region,  and  d)  the  RC  time  constant  of 
the  device  combined  with  external  load  resistance.  In  the 
following  we  will  discuss  these  mechanisms  in  detail.  The 
diffusion  time  can  be  minimized  if  a reach-through 
avalanche  photodiode  (RAPD)  structure  is  employed  [33]  as 
shown  in  Fig.  2.12.  In  this  case  the  response  time  of  the 
internal  APD  device  is  limited  by  only  three  mechanisms: 

a)  the  time  for  carriers  to  complete  the  process  of 
multiplication,  t&; 

b)  the  electron  transit  time  across  the  drift  region 
(Wd) / tdfe  = Wd/vse/  where  vse  represents  the  electron 
saturation  velocity; 

c)  the  transit  time  of  the  last  holes  across  the 
avalanche  region  (Wa)  and  the  drift  region, 

_ Wa  + Wd 
^dh  - ,,  , 

vsh 


where  vSh  denotes  the  hole  saturation  velocity. 

The  avalanche  build-up  time,  tA,  is  a function  of  the 
ratio  of  the  ionization  coefficients  k and  and  given  by 


[34] 


tA  — Mk 


wA 


vSe 


0 < k < 1 


(2.160) 
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n 


n+ 

Substrate 


(a) 


Figure  2.12  Reach-through  avalanche  photodiode 

(a)  Device  configuration 

(b)  Doping  profile 

(c)  Field  distribution 

(d)  Absorption  of  incident  light. 


64 


The  overall  response  time  becomes 

_ Wd  + MkW/y  Wd  + Wa 
vse  vsh 


(2.161) 


The  bandwidth  is  characterized  by  the  3dB  cutoff  frequency 
and  is  approximately  given  by  [34] 


^ 0.44 

f(-3dB)  = — ~ 


(2.162) 


where  0.44  is  a numerical  factor. 

Equations  (2.160)  and  (2.161)  show  that  increasing 
the  multiplication  factor  results  in  decreasing  the 
bandwidth.  That  is,  an  APD  exhibits  a limited  gain- 
bandwidth  product . 

The  effect  described  by  Eq.  (2.161)  gives  another 
reason  for  seeking  materials  with  k « 1 for  devices  in 
which  the  electrons  initiate  the  avalanche. 

2. 3. 1.4.  Quantum  Efficiency 

The  quantum  efficiency  of  the  p-i-n  structure  is 
defined  by  the  number  of  electron-hole  pairs  generated  per 
incident  photon  [35] : 


^ = iph/q 
Popt/fiv 


(1  - r)  (1 


p-aw 

— ) 

1 + aLp 


(2.163) 


where 

Ip  = photogenerated  current  by  the  absorption  of 
light . 
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Popt  = incident  optical  power 
v = light  frequency 
r = reflection  coefficient 
a = absorption  coefficient 
W = depletion-layer  width 

Lp  = minority  carrier  (hole)  diffusion  length 
Eq.  (2.163)  shows  that  one  of  the  key  factors  that 
determines  the  quantum  efficiency  is  the  absorption 
coefficient.  Fig.  2.13  shows  the  measured  intrinsic 
absorption  coefficients  a for  various  photodetector 
materials  [36]  . Since  a is  a strong  function  of  the 
wavelength,  the  useful  wavelength  range  for  detection  is 
limited  for  a given  semiconductor  photodiode.  The  long 
wavelength  cutoff  Xc  is  established  by  the  energy  gap  of 
the  semiconductor.  For  example, 

i 1.24 

/'c  Eg (=  0.67)  1,8  M*11  for  Ge 

1 . 24 

K:  = y:'i4  ~ 1 . 1 |im  for  Si 
1 . 24 

= l 43  ~ 0 . 9 Jim  for  GaAs 

For  wavelengths  larger  than  A.c,  the  values  of  a are  too 
small  to  give  appreciable  absorption.  The  short  wavelength 
cutoff  of  the  photodiodes  comes  about  because  the  values 
of  a for  short  wavelengths  are  very  large  (>  105cm_1) , and 
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Figure  2.13  Optical  absoption  coefficients  for 
various  photodetector  materials. 
(Ref.  36) . 
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the  radiation  is  absorbed  very  near  the  surface  where  the 
recombination  is  fast.  Thus  the  photocarriers  recombine 
before  they  are  collected  in  the  p-n  junction. 

Another  result  we  can  deduce  from  Eq.  (2.163)  is  that 
to  obtain  a high  quantum  efficiency,  the  depletion  layer 
must  be  wide  enough  to  ensure  that  aw  » 1.  However,  the 
transit  time  delay  may  be  considerable.  Practically,  the 
depletion  width  is  chosen  equal  to  about  2/a.  At  this 
depth  the  incident  optical  power  falls  to  a fraction  1/e2 
of  its  value. 

The  second  key  factor  in  the  expression  for  the 
quantum  efficiency  is  the  reflection  coefficient  r.  In 
most  of  the  materials  suitable  for  photodiodes,  the 
discontinuity  in  the  refractive  index  going  from  air  (n  = 
1)  to  the  semiconductor  (n  ~ 3.5)  leads  to  a reflection 
coefficient  r of  0.3  at  normal  incidence.  This  reflection 
loss  can  be  eliminated  by  depositing  an  appropriate 
transparent  antireflection  film  coating  on  the  photodiode 
surface  [37] . Silicon  nitride,  with  a refractive  index  in 
the  range  n = 1.8  ~ 2.0,  can  be  used  for  this  purpose.  The 
film  thickness  must  be  equal  to  a quarter  wavelength  of 
the  light  in  the  dielectric. 

For  practical  purposes,  instead  of  the  quantum 
efficiency,  the  reponsivity  R is  used  and  it  is  defined  by 
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r = TPh  = sn 

popt  hv 


(2.164) 


which  is  the  ratio  of  photocurrent  Iph  to  incident  optical 
power  P0pt- 


2. 3. 1.5.  Dark  Current 

The  reduction  of  the  dark  current  is  important  since 
it  limits  the  minimum  detectable  power.  In  an  APD  the  dark 
current  has,  in  general,  three  components:  a)  the 
generation  and  recombination  current,  b)  the  minority 
diffusion  current,  and  c)  the  tunneling  current  in  the 
avalanche  region. 

For  low  temperatures  the  dark  current  is  governed  by 
generation-recombination  currents  which  are  proportional 
to  the  depletion  width,  which  in  turn  depends  on  the  bias 
voltage.  However,  in  the  high  temperature  range,  above 
room  temperature,  the  dark  current  is  dominated  by 
diffusion  currents. 

In  another  classification  the  dark  current  can  be 
divided  into  two  components:  a)  dark  current  that  is 
multiplied,  Idm  by  the  avalanche  effect;  and  b)  dark 
current  that  is  not  multiplied,  IdO/  and  which  flows 
through  the  periphery  of  the  guard  ring.  In  terms  of  shot 
noise  we  get  [38] 

= 2qI3{M2  (Iph  + Idm)F  + Ido} 


(2.165) 
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where  B is  the  bandwidth  and  Iph  is  the  photosignal 
current.  If  M is  large,  Idc>  can  be  neglected  and  the  dark 
current  is  governed  by  Idm  and  a linear  relationship  is 
obtained  between  the  dark  current  and  the  multiplication 
factor  (Fig.  2.14)  . Thus,  Idm  is  determined  by 
extrapolating  the  linear  relationship  to  unity 
multiplication . 

In  the  following  we  will  discuss  the  three  components 
of  the  dark  current  in  more  detail . 

a)  The  Generation  and  Recombination  Current.  Tcr-r 

The  value  of  Ig-r  is  proportional  to  the  volume  of 
the  depletion  region  and  inversely  proportional  to  the 
effective  carrier  lifetime.  Ig_r  is  given  by  [37] 

T qn-jAW  qv 

Ig-r  = — — — [1  - exp(-^jr)]  (2.166) 

where  ni  is  the  intrinsic  carrier  concentration.  The 
magnitude  of  this  current  is  strongly  affected  via  Te  by 
crystal  imperfections,  impurities,  etc.,  within  the 
depleted  region.  High  quality  of  defect-free  material  is 
required  to  reduce  Ig_r [37]  . 

b)  The  Diffusion  Current 

The  diffusion  current  Idif  arises  from  thermally 
generated  minority  carriers  diffusing  into  the  depletion 
region  from  the  surrounding  undepleted  p and  n regions. 
The  contributions  of  diffusion  current  depend  on  the 


Dark  Current  ( jiA) 
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Multiplication  Factor  (M) 


Figure  2.14  A plot  of  the  total  dark  current 

versus  multiplication  factor.  Note 
that  separation  of  IDM  is  obtained 
from  the  graphical  method. 
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details  of  the  photodiode  geometry  and  can  be  calculated 
from  [37] . 


Is  = 


<^> 1/2 


p region 


(2.167) 


Is  = qn?  (^)1/2(^-)  : n region  (2.168) 

where 

Dn  = n-type  minority-carrier  diffusion  constant 
Dp  = p-type  minority-carrier  diffusion  constant 
xn  = n-type  minority-carrier  lifetime 
Xp  = p-type  minority-carrier  lifetime 
Ap  = area  of  depletion  region  bounding  the  p-region 
An  = area  of  depletion  region  bounding  the  n-region 
Na  = acceptor  concentration 
Nd  = donor  concentration. 

Because  the  diffusion  current  contribution  to  the  dark 
current  is  proportional  to  n?  and  the  Ig-r  contribution  is 

proportional  to  ni,  Idif  becomes  dominant  in  materials 
with  small  bandgaps  and  at  high  temperatures, 
c)  The  Tunneling  Current  At  sufficient  high  fields, 
tunneling  can  take  place  as  a band-to-band  process  or  via 
traps  or  defect  levels  in  the  forbidden  energy  gap.  The 
expression  for  band-to-band  tunneling  in  a uniformly  doped 
depletion  region  may  be  written  as  follows  [37] 


Itun 


y-A-exp  [- 


*&/2e3/2 

qhEm 


(2.169) 
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where  mQ  is  the  free  electron  mass  and  Em  is  the  maximum 
junction  field.  The  parameter  0 = a(^p)1//2'  where  me  is 

the  electron  effective  mass  and  a ~ 1,  and  depends  on  the 
detailed  shape  of  the  tunneling  barrier.  The  factor  y is 
different  for  tunneling  via  traps  in  the  bandgap  and  for 
band-to-band  tunneling.  In  the  latter  case. 


v _ , 2nu  qEmV 

Y \\,2)  47l2h2 


(2.170) 


From  the  exponential  term  in  Eq.  (2.169),  it  can  be 
seen  that  the  tunneling  current  increases  with  temperature 
because  of  the  decrease  of  Eg  with  increasing  temperature. 
The  tunneling  current  near  breakdown  becomes  a serious 
problem,  especially  for  long  wavelength  APD  (where  Eg  is 
small) . To  reduce  the  parasitic  effect  of  the  tunneling 
current  near  breakdown  voltage,  an  InGaAs/InP  hetero- 
junction structure  has  been  proposed  [33,40]  as  we  will 
explain  in  the  next  section. 


2.3. 1.6.  Avalanche  Photodiode  Structure 

In  optical  fiber  communication,  a detector  is 
required  to  have  (i)  a low-noise  performance,  (ii)  a low- 
operation  voltage,  (iii)  a stable  gain,  and  (iv)  high- 
speed performance.  The  reach-through  structure  is  designed 
to  meet  all  these  requirements  [33] . 
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The  reach-through  avalanche  photodiode  ( RAPD ) 
structure  shown  in  Fig.  2.12  has  the  following  advantages. 

1)  The  absorption  in  the  p+  region  is  small  and  the 
slow  diffusion  of  carriers  to  the  high-field  region  is 
minimized. 

2)  The  main  absorption  occurs  in  the  high-field  thick 
n-region  where  a high  quantum  efficiency  and  a small 
transit  time  can  be  achieved. 

3)  The  p-n+  junction  is  guarded  by  the  surrounding  n 
region  so  that  uniform  avalanche  gain  can  be  achieved  in  a 
uniform  junction  without  edge  breakdown. 

4)  The  doping  of  the  p-n+  junction  where  avalanche 
multiplication  occurs  can  be  optimized  to  obtain  the 
lowest  possible  noise  and  to  obtain  a slowly  varying 
multiplication  factor  with  bias  voltage. 

In  III -V  compound  materials  (direct  bandgap 
semiconductors),  a "soft"  breakdown  of  the  dark  current 
due  to  tunneling  is  typical  [41] . To  solve  this  problem 
Nishida  [42]  made  use  of  the  heterostructure  shown  in  Fig. 
2.15.  In  this  structure  the  GalnAs  layer  forms  a narrow 
bandgap  region  designed  to  absorb  the  incident  photons. 
The  InP  layer  is  a separate  wide  bandgap  region,  included 
to  support  avalanche  multiplication.  Any  APD  structure 
having  these  separate  regions  has  become  known  as  a 
separate  absorption  and  multiplication  APD  (SAM-APD) . 


I 
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p+-  InP  inP  n-  InGaAs 


(a) 


Figure  2.15  Energy  band  diagrams  of  an  InGaAs-InP 

photodiode  with  separate  absorption  and 
multiplication  regions  (Ref.  32) 

(a)  zero  bias 

(b)  reverse  bias. 
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The  narrow  bandgap  region  must  be  sufficiently  thick 
to  absorb  all  the  incident  photons  in  order  to  yield  a 
high  quantum  efficiency  and  must  have  a sufficiently  low 
doping  concentration  for  it  to  be  fully  depleted  with  a 
modest  applied  field.  The  thickness  and  doping  level  of 
the  wide  gap  material  are  adjusted  so  that  the  electric 
field  in  the  narrow  bandgap  material  never  exceeds  the 
value  which  causes  significant  tunneling  current. 

The  drawback  of  this  structure  is  that  this  APD  has  a 
poor  time  response  due  to  a pile-up  of  holes  at  the 
discontinuity  in  the  valence  band  at  the  InGaAs/InP 
interface.  Capasso  [43]  overcame  this  drawback  by  using 
superlattice  layers  at  the  interface  grading  out  the 
discontinuity. 

2.3.2.  Earameter  Measurements 

The  first  measurements  in  which  the  ionization 
coefficients  a and  (3  were  separately  determined  from 
photomultiplication  measurements  on  a P+N  junction  were 
reported  by  Wul  and  Shotov  [44].  Their  analysis  assumed  a 
constant  value  of  the  ratio  k = (3/a. 

Lee  et  al . [45]  showed  that  if  the  measurements  of 
the  electron  and  hole  multiplication,  Mn  and  Mp,  are  made 
in  the  same  junction,  then  determination  of  a and  |3  is 
possible  from  the  relationship 
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w 

f M 

J (a-p)dx  = Infer)  (2.171) 

0 * 

To  determine  accurately  Mn  and  Mp,  several  experimental 
conditions  must  be  met  [46] : 

1)  Pure  electron  and  pure  hole  injection  must  be 
obtained  in  the  same  junction.  This  is  best  achieved  by 
using  photoexcitation. 

2)  The  photocurrent  without  avalanche  gain  must  be 
measured  precisely. 

3)  The  electric  field  variation  in  the  junction 
should  be  sufficiently  gradual  so  that  a and  P can  be 
assumed  to  be  only  functions  of  electric  field. 

4)  The  electric  field  profile  must  be  accurately 
known . 

5)  The  avalanche  gain  must  be  uniform  across  the 
active  area  of  the  device. 

There  are  several  device  structures  that  allow  for 
the  accurate  determination  of  oc  and  p.  Here  we  introduce 
the  p-i-n  diode  structure  which  is  the  simplest  and 
satisfies  all  of  the  requirements  just  listed. 

Fig.  2.16  shows  the  p-i-n  diode  configuration  and  its 
electric  field  distribution.  Pure  electron  and  hole 
injection  can  be  obtained  because  the  incoming  light  is 
completely  absorbed  in  the  p+  or  n+  region.  Since  the 
electric  field  is  constant  in  this  structure  (E  = v/W, 
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Figure  2.16  A p-i-n  photodiode  and  its 
electric  field  profile. 
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where  V is  the  total  voltage  across  the  intrinsic  region) , 
the  ionization  coefficients  are  constant  and  given  by 


1 r Mn  - 1 , , ,Mni 
a = » [ — I"—  ] In  (jr-) 


W L Mn  - Mp 


(2.172) 


P = 


_ 1 r MP  - 1 


.Mr 


W t Mp  - Mn  ^ ln(Mn) 


(2.173) 


In  the  following  we  will  derive  equations  (2.172)  and 
(2.173)  . 

From  Eqs.  (2.153)  and  (2.154)  we  find 


w 


Mn 


exp[  J (a-p)dx] 
0 


w 


w 


w 


w 


x {1  - J p exp[  J (a-P)dx']dx} 


1 - J P exp[  J (a-P)dx']dx 


X 


X 


w 


= exp[  J (a~P)dx] . 


(2.174) 


Therefore, 


W 


f M 

0(w)  = J (a-p)dx]  = lnfcjj1) 

0 ^ 


(2.175) 


which  is  the  same  as  Eq.  (2.171).  Eq.  (2.174)  can  be 
rewritten  as 


a 1 ,Mn. 

“ - P = w ln(5^> 


(2.176) 
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From  Eq.  (2.154) 


w 


w 


Mp  - 1 = 


exp [-  J (cc-p)dx]  - 1 + J a exp [-  J (a-P)dx']dx 

o o o 

W X 

1 - J aexp[-J  (a-J3)dx']dx 
o o 

(2.177) 

By  subtracting  Eq.  (2.153)  from  Eq.  (2.154),  we  have 


w 


Mp  - Mn  = 


exp[-J  (a-P)dx]  - 1 

o 

W X 

1 - J a exp [-  J (a-P)dx']dx 
0 0 


(2.178) 


The  denominators  of  Eqs.  (2.177)  and  (2.178)  are  the  same, 
Thus 


w 


w 


x 


Mp  ~ 1 

Mp  - Mn 


exp[-J  (a-P)dx]  - 1 + J aexp[-J  (cx-P)dx']dx 
0 0 0 


w 


= 1 + 


= 1 - 


exp[-  | (a-p)dx] 

0 

W X 

J a exp[-J  (a-p)dx] 

o g 

W X 

- J (a-P)exp[-J  (a-P)dx']dx 

o o 

a 


(a  - p) 

(2.179) 

By  rearranging  Eq.  (2.179)  in  terms  of  p,  we  obtain 


__ 

P - a 
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P *P  - Mn  ‘ 


Substitution  of  Eq.  (2.176)  yields 


P = W [ M^-  M-  ] ln  ' 


‘n 


M 


ln 


(2.180) 


(2.181) 


In  a similar  way,  we  have  the  expression  for  a as 


« = w [ J?n.  i ] ln^T)  Q.E.D. 


W L Mn  - Mr 


Fig.  2.17(a)  shows  two  practical  device  structures 
that  allow  for  very  pure  optical  injection  of  electrons 
and  holes.  These  have  become  the  standard  device 
configurations  for  ionization  coefficient  measurements 
[37]  . 


In  the  first  structure  note  that  a window  is  etched 
in  the  substrate  to  within  a few  tens  of  micrometers  from 
the  junction  in  order  to  obtain  pure  electron  injection 
using  strongly  absorbed  light.  However,  in  materials  with 
a short  minority-carrier  diffusion  length,  many  of  the 
electrons  optically  excited  in  the  region  of  the  substrate 
just  above  the  etched  window  recombine  before  being 
collected  by  the  junction.  The  near-bandgap  radiation 
generated  by  this  recombination  can  be  reabsorbed  by  the 
Franz-Keldysh  effect  in  the  space-charge  region.  This 
reabsorption  produces  a contamination  of  the  electron- 
injected  current,  and  if  the  width  of  the  avalanche  region 
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is  appreciable,  pure  electron  injection  is  impossible. 
This  effect  was  first  observed  by  Bulman  et  al.  [47]. 

Capasso  [48]  introduced  a modification  of  the  basic 
structure  (Fig.  2.17(b))  to  eliminate  this  source  of 
contamination.  A thin  (<  1 jim)  buffer  layer  (p+-InGaAsP) , 
composition  lattice  matched  to  the  substrate,  is  inserted 
between  the  n-layer  and  the  substrate. 

The  window  is  etched  with  a selective  etchant  that 
stops  at  the  buffer  layer.  This  thin  layer  with  a smaller 
bandgap  completely  absorbs  the  incident  light;  yet  the 
thickness  is  substantially  smaller  than  the  minority 
carrier  diffusion  length,  so  that  the  amount  of 
recombination  is  negligible. 

Fig.  2.18  shows  an  experimental  set-up  for  the 
measurements  of  the  ionization  rates  [37] . The 
photocurrent  is  measured  via  phase-sensitive  detection.  An 
optical  fiber  is  used  to  couple  the  light  source  to  the 
photodiode.  This  minimizes  any  contamination  of  the 
electron  photocurrent  from  stray  light  incident  on  the 
edge  of  the  mesa. 
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Figure  2.17  Configuration  of  practical  device 
structures  (Ref.  6) 

(a)  a simple  structure  used  for  the 
parameter  measurements. 

(b)  an  improved  structure  for  pure 
electron  or  hole  inection. 
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Figure  2.18  Typical  experimental  set-up  used  in 

photomultiplication  measurements  (Ref.  6). 


CHAPTER  III 

TWO  SIMPLE  PROOFS  OF  CAPASSO'S  EXECESS  NOISE  FACTOT  FN 
OF  AN  IDEAL  N-STAGE  STAIRCASE  MULTIPLIER 

We  begin  with  a simple  proof  of  Capasso's  excess 
noise  factor  [49]  of  an  ideal  N-stage  staircase  multiplier 
by  applying  the  variance  theorem  [50]  to  an  ideal  one- 
stage-staircase  multiplier.  The  applicability  of  the 
variance  theorem  to  a branching  process  is  not  immediately 
obvious  since  correlations  between  the  various  generations 
occur.  Yet  the  overall  result  indicates  the  validity  of 
the  variance  theorem  as  was  shown  explicitly  by  Van  Vliet 
[51] . "Ideal"  means  that  only  the  electrons  contribute  to 
the  device  operation.  The  noise  process  of  such  a one- 
stage  amplifier  can  be  described  in  the  following  way.  We 
write 

ni 

mi  = £xk  (3.1) 

k=l 

where  mi  and  ni  are  the  fluctuating  output  and  input 
particle  numbers,  respectively,  and  Xk  is  a stochastic 
event,  in  this  case  the  impact  ionization  process  at  the 
staircase  interface.  Then,  Xk  = 1 (no  impact  ionization) 
with  probability  (1  - A,)  and  Xk  = 2 with  probability  X. 
Consequently 
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Xk  = X = Mi  = 1 + A, 


(3.2a) 


and 


var  Xk  = var  X = (1  - A,)  -A 


(3.2b) 


In  Eq.  (3.2a),  Mi  is  the  average  gain  of  a single  ideal 
stage.  By  applying  the  variance  theorem  to  (3.1),  we 
obtain 

o 

var  mi  — X var  ni  + ni  var  X (3.3a) 

Since  we  assume  a Poissonian  distribution  for  n i, 
Eq. (3.3a)  reduces  to 

var  mi  = ni(X2  + va  X}  (3.3b) 


The  excess  noise  factor  Fi  of  a one-stage  multiplier 
is  defined  by 


Fi 


noise  in  output  particle  number 

amplified  input  noise  if  Xk  did  not  fluctuate 


var  mi 
X var  ni 


(3.4a) 


Upon  substituting  Eqs.  (3.2a)  and  (3.3b)  into  Eq.  (3.4a), 
we  find  for  Fi 


X2  + var  X 
X2 


„ var  X 
1 + 

X2 


1 + 


A(l-A) 

(1+A)2 


Fi  = 


(3.4b) 
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If  Mn  is  the  noise  measure  of  the  multiplier  [52], 
defined  as 

(1  - X) 


Fj  - 1 


*n  — 


1 - 1/Mi  (1  + ^ 


(3.5) 


then  the  excess  noise  factor  Foo  of  an  infinite  number  of 
identical  stages  in  cascade  is  (idln  + 1)  . If  the  infinite 
number  of  stages  is  divided  into  groups  of  N subsequent 
stages,  each  with  an  excess  noise  factor  Fn  and  a gain  Mn 
= (1  + X) N,  then  Eq.  (3.5)  may  also  be  written  as 

(1  ~ X)  Fn  ~ 1 


Mn  = 


(l  + X) 


1 - 1/Mn 


(3.6) 


Consequently 


FN  = 1 + 


1 - X\ 

1 + Xj 


[1 


(1  + A.)-N] 


(3.7) 


which  is  Capasso's  result  [49], 

The  second,  even  simpler  proof  uses  Friiss's  formula 
[53]  for  the  excess  noise  factor  Fn  of  N identical  single 
stages,  each  of  excess  noise  factor  Fi  and  gain  Mi. 
According  to  Friiss 

Fi  - 1 Fi  - 1 

FN  = 1 + (Fi  - 1)  + — “ + • • • + — 

Mn  (Mn)N-1 

= 1 + — ~ ~ 1 [1  - (Mn)  -»]  (3.8) 

1 - (Mn)"1 

Substituting  for  Fi  and  Mi  yields  again  Eq. (3.7). 
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The  introduction  of  the  noise  measure  j$ln  or  the 
application  of  Friiss's  formula  Eq.  (3.8)  puts  the 
calculations  of  individual  situations  on  a general 
schematic  basis  that  is  easy  to  apply  to  this  and  to  other 
systems;  for  example,  Eq.  (3.8)  applies  also  to  secondary 
emission  multipliers  and  it  connects  signal  amplifier 
theory  with  particle  multiplier  theory. 

Eq.  (3.7)  agrees  with  the  results  presented  in  Van 
Vliet's  single-carrier  multiplication  paper  [51].  However, 
in  this  paper  no  use  is  made  of  the  variance  theorem,  but 
rather  the  branching  problem  is  solved  by  an  iteration 
process  that  involves  the  consideration  of  correlation 
between  the  branches  of  offspring.  Teich  and  Matsuo  [54] 
have  shown  that  Van  Vliet's  method  and  Capasso's  method 
lead  to  the  same  results.  When  the  holes  also  contribute. 
Van  Vliet's  two  carrier  multiplication  process  [27]  must 
be  applied. 


CHAPTER  IV 

NONPARABOLIC I TY  AS  A TOOL  IN  BAND-GAP  ENGINEERING 

4.1.  Introduction 

The  study  of  impact  ionization  thresholds  in 
semiconductors  was  shown  to  be  equivalent  to  the  study  of 
the  "most  probable'  Auger  transitions  in  the  same  system 
some  time  ago  [55] , and  the  latter  study  goes  back  even 
further  [56] . Many  papers  on  threshold  calculations  have 
been  published  since  that  time.  With  increasing  interest 
in  low-noise  avalanche  photodiodes  and  impact  ionization 
across  energy-band  discontinuities  in  heterostructures 
using  III-V  semiconductor  compounds  and  related  materials, 
it  is  of  interest  to  reexamine  what  is  known  about 
threshold  calculations  with  special  reference  to  band 
nonparabolicity . For  processes  involving  up  to  three 
parabolic  bands  standard  formulas  are  of  course  available 
[57]  . We  note  early  calculations  of  thresholds  for 
nonparabolic  bands  both  analytically  [57]  and  numerically 
[58,59].  Graphical  constructions  have  also  been  used  [60- 
62].  For  later  work,  recent  reviews  [62,6]  may  be 
consulted. 
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4.2.  Derivation  of  a Formula  for  the  Impact - 

IonizationThreshold  Energy  with  Nonparabolic 
Conduction  and  One  or  Two  Parabolic  Valence  Band 

As  we  found  in  chapter  II,  momentum  conservation 

principle  yields 


where  ki  and  ki  represent  wave  vectors  of  initiating 

particle  and  offspring  particle  1,  respectively.  If  we  use 
the  definition  of  |i  = this  equation  reduces  to 


We  assume  that  an  initiating  electron  is  located  in 
the  conduction  band  of  the  simple  Kane  type.  Then 


where  me  is  the  effective  mass  at  the  conduction-band 
minimum  and  a'  is  a nonparabolicity  parameter  for  a simple 
Kane  band.  For  the  threshold  condition  Eq.  (4.2)  can  be 
rewritten  as 


Since  this  is  a simple  quadratic  equation,  we  solve  for  Et 


(2.19) 


(4.1) 


(4.2) 


(4.3) 


and  obtain 
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Et  ' - to  + 2^  + 4“’A 


(4.4) 


where 


A 


2 2 
ft  kj 

2 me 


(4.5) 


Alternatively,  threshold  energy  can  be  obtained  by  energy 
conservation  principle,  which  is 


ft2ki  ft2k2 
Et  * — + — 


tl2kj 

- <-  > + 


(4.6) 


where  we  assumed  that  two  conduction-band  states  (ki  and 
k2)  and  one  valence-band  state  (k3)  are  treated  as 
parabolic.  This  is  possible  as  they  are  states  of  low 
energy.  It  is  this  physical  approximation  which  makes  our 
problem  tractable  by  analytical  methods.  This  is  not 
possible  if  one  treats  all  three  conduction-band  states  as 
nonparabolic  [63] . 

Upon  substituting  Eq.  (4.1)  into  (4.6),  we  obtain 


Et 


,22 
ft  kj 

2 i% 


[- 


2 2 
ftki 

2 me 

1 + 2|i 


(1  +2|i)  2 

r \L 


+ 


1 + 2 \L 
+ Eq 


(1  +2(1) 
']  + Eq 


2]  + Eg 


(4.7) 


By  equating  Eq.  (4.4)  and  (4.7)  with  some  rearrangements, 
we  have 


91 


A2  (f^) 2 + A[8a,2Eb  4a'  ] + 4a,24  + 4a'EG  = 0 

If  we  use  substitutions  such  as 

1+2*1  S 

then,  Eq.  (4.8) 


U±tL=  e 

l-2|i  ' 


reduces  to 


(4.8) 


A2  (4a' 2<{)2)  + A [8a' 2^  - 4a' 6]  + (4a'2Eg  + 4(x'Eg)  =0  (4.9) 


which  is  a quadratic  equation  form  for  A.  If  we  solve  this 
equation  for  A,  we  obtain 


A = 


l-4a'EG 


M-  (1+2*1) 
(1+M-) 2 


ex' Eg 

♦ 


(4.10) 


Finally, 

have 


we  substitute  Eq.  (4.10)  into  (4.7).  Then,  we 


where  it  turns  out  that  only  the  bottom  sign  in  the 
bracket  has  physical  meaning. 

Ridley's  formula  (4.6)  can  be  easily  derived  if  we 
use  the  condition  such  that 
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(4.12) 


This  condition  yields  an  approximation  of  the  square  root 
term  between  the  brackets  of  Eq.  (4.11).  That  is. 


Using  this  first  order  approximation,  Ridley's  formula  is 
obtained  as 


If  one  is  concerned  with  GaAs  and  uses  p.  ~ 0.067/0.45 
= 0.149,  (XEq  ~ 1,  Eq  ~ 1.424  eV,  this  formula  gives  Et  ~ 
1.84  eV,  which  is  a distinctive  improvement  on  the  value 


which  one  would  find  for  parabolic  bands,  if  one  notes 
that  the  accepted  value  lies  near  2 eV  [6]  . Using  the 
values  for  GaAs,  our  new  formula  Eq.  (4.11)  yields  the 
result  Et  = 1.95  eV,  in  excellent  agreement  with  the 
accepted  value. 

Note  that  Eqs.  (4.11)  and  (4.14)  share  as  common 
characteristic  that  they  reduce  to  Eq.  (4.15)  in  the  limit 
OtEG  — » 0,  as  is  of  course  required. 


(4.14) 


[56,  57] 


l + 2ll 

Et  = "l  + (1  eG  ~ 1 • 61  eV 


(4.15) 
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4.3.  Further  Discussions  Concerning  the 
Nonparabolicity 

Formula  (4.11)  is  compared  with  (4.14)  in  Fig.  4.1, 
showing  that  the  new  formula  yields  higher  threshold 
energies.  There  is  a cutoff  when  the  square-root  becomes 
imaginary  and  formula  (4.11)  suggests  that  impact 
ionization  ceases.  However,  for  these  higher  Et  values  the 
approximation  just  discussed  begins  to  fail  as  the  two 
low-energy  states  1 and  2 (Fig.  4.2(b))  move  higher  into 
the  conduction  band  and  begin  to  feel  the  nonparabolicity. 
Thus  the  cutoff  value  [L  = (ic  suggested  by  Eq.  (4.5)  which 
occurs  when 


a' EG 


(1  + M-)2 
4*1(1  + 2 n) 


(4.16) 


or 

V (a'EG)2  + a'Es  - (cc'Eg  - y) 

Mc  = 2 (“'Eg  > 0 > (4.17) 

can  be  misleading.  To  establish  a more  reliable 
expression,  we  look  for  the  condition  that  the  line  on  the 
E-k  diagram  representing  the  downward  transition  in  the 
conduction  band  cannot  intersect  the  valence  band.  This 
means  that  the  valence-band  asymptotes  must  be  steeper 
that  the  conduction-band  asymptotes  (Fig.  4.2).  The 
condition  for  this  turns  out  to  be  simply 


ErfeV)  ErteV) 
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CxEo 

(a) 


(b) 


Figure  4.1  Threshold  energy  as  a function  of 

(a)  aEg  (using  }i  = 0.149)  and  (b)  |i  (using 
aEg  = 1.0)  for  GaAs . The  present  result  is 
represented  by  a solid  line,  dashed  line 
represents  Ridly's  result,  and  the 
experimental  value  is  shown  as  a circle. 
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Figure  4.2  Quadruplets  of  states  at  the  trheshold  of 
impact  ionization  at  given  aEg  = 1.  The 
effect  of  increasing  m is  shown.  The 

horizontal  axis  gives  ak/p  where  a =5.65&. 
(a)  x>  = 0.2  nc,  (b)  (i  = \iCf  (c)  [i  = 1.2  |i*. 
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M-  < M- 


(4.18) 


Thus  for  a' Eq  ~ 1,  impact  ionization  is  still  possible  in 
the  range  (see  Fig.  4.3) 


It  ceases  for  |i  > |i*,  provided  the  Kane  band  model  still 
holds . 

Thus  the  new  formulas  (4.11)  and  (4.18)  show  that  one 
has  a way  of  "switching  off"  this  particular  impact 
ionization  process,  given  only  a band  of  the  form  Eq. 
(4.2)  . In  a band-gap  engineering  context  if  one  needs 
strong  impact  ionization  one  is  looking  for  small  (ivalue- 
this  is  already  suggested  by  the  older  formula  (4.15).  One 
can  on  the  other  hand  decrease  it,  or  indeed  switch  it  off 
altogether,  by  looking  for  large  enough  values  of  [i.  This 
latter  feature  is  not  even  hinted  at  by  the  earlier 
formulae  (4.14)  or  (4.15).  In  view  of  the  occurrence  of 
more  complicated  band  structures  (see  for  example  Ref. 
12)  , the  problem  of  identifying  other  structures  also 
displaying  switch  off  is  of  interest  and  is  under 


|Ic  < \l  < |I 


i.e.,  0.26  <|l<  0.5 


consideration . 
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Figure  4.3  Electron/hole  effective  ma; 

M-c  and  |i*as  a function  of  ( 


ratios 

g* 


Chapter  V 

THE  THEORY  OF  TRAPPING  AND  DETRAPPING  OF  ELECTRONS 
IN  MULTIQUANTUM-WELL  (MQW)  HETEROSTRUCTURES 

5.1.  Introduction 

It  is  well  known  that  noise  is  associated  with  the 
randomness  of  the  charge  transport  mechanism.  If  there  is 
a fluctuation  inside  an  electronic  device,  the  electric 
signal  at  the  external  terminal  will  fluctuate  as  a 
result.  One  of  the  fluctuation  mechanisms,  the  carrier 
generation  and  recombination  (g-r)  process,  is  quite  a 
common  source  of  carrier  density  fluctuations  in 
semiconductor  devices.  This  g-r  process  takes  place  via 
trap  energy  levels  located  in  the  forbidden  band-gap. 

The  noise  associated  with  the  random  filling  and 
emptying  of  the  traps  is  called  trapping  noise  and  has 
been  studied  for  many  decades.  Gisolf  [65]  used  a Fourier 
Transform  method  to  find  the  power  spectrum  of  a current 
pulse  resulting  from  carrier  detrapping  and  whose  duration 
is  determined  by  the  time  between  the  release  of  a carrier 
from  a trap  and  the  subsequent  capture  by  a trap.  He 
assumed  each  event  to  be  independent  and  a fixed  carrier 
lifetime.  The  latter  assumption  is  not  quite  right  because 
the  carrier  lifetime  in  reality  is  distributed.  Therefore, 
Bernamont  [66]  and  van  der  Ziel  [67]  took  into  account  a 
probability  function  for  the  lifetime  of  free  electrons. 
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assumed  it  to  be  Poissonian,  and  came  up  with  a Lorentzian 
type  spectrum.  These  researchers  considered  only  the 
lifetime  of  the  free  carriers.  Machlup  [68]  however,  took 
also  into  account  the  mean  time  carriers  are  bound  in 
traps.  Thus  he  compared  the  trapping  and  detrapping  with  a 
random  telegraph  signal  which  switches  between  two  states 
randomly.  If  the  mean  free  time  is  relatively  small 
compared  to  the  mean  bound  time,  his  result  approaches 
Bernamont's  result. 

More  general  approaches  have  been  developed  later  by 
Van  Vliet  et  al.  [69,70].  They  were  interested  in  how  to 
incorporate  all  the  information  concerning  the  microscopic 
and  macroscopic  electron-trap  behavior  into  the  basic 
equations  which  govern  device  operation.  This  so-called 
collective  approach  consists  basically  of  two  parts. 
First,  the  internal  sources  of  the  fluctuations  due  to  g-r 
processes  and  diffusion  are  accounted  for.  Secondly,  the 
response  function  which  relates  the  external  terminal 
noise  to  the  internal  distributed  sources  should  be 
calculated.  These  two  parts  are  already  linked  together  in 
the  Fourier  Transformed  transport  equations.  An  elegant 
formalism  was  presented  by  Van  Vliet  et . al  and  named  the 
Transfer  Impedance  (T.I.)  method.  In  general,  an  analytic 
solution  of  the  differential  transport  equations  is 
possible,  though  in  many  cases  it  becomes  very 


100 


complicated.  We  follow  the  collective  T.I.  approach  to 
solve  our  problem. 

G-r  noise  spectra  obtained  by  solving  the  transport 
equations  can  be  characterized  by  basically  two 
parameters,  namely  by  the  low-frequency  plateau  values  of 
the  spectra  and  by  the  cutoff  frequencies  of  the 
Lorentzians.  Both  of  the  parameters  depend  on  the 
activation  energy,  density  and  capture  coefficient  of  the 
traps . 

From  our  noise  measurements,  we  observed  substantial 
sub-shot  noise  in  the  high-frequency  regime  of  some 
devices,  indicating  that  carrier  flow  in  the  MQW 
structures  is  severely  modified  by  the  trapping  and 
detrapping  in  the  quantum-wells. 

To  understand  the  sub-shot  noise  associated  with  the 
trapping,  we  look  at  the  carrier  flow  in  the  presence  of 
traps.  Under  reverse  bias  condition,  carrier  motion  is 
limited  by  the  transit  time  across  the  depletion  region. 
If  traps  are  present,  carrier  arrival  is  delayed  because 
carriers  spend  some  time  in  the  traps  until  they  are 
released  to  resume  their  flight.  Since  current  is  charge 
flow  per  unit  time,  trapping  results  in  a current  decrease 
by  a factor  equal  to  a time  constant  ratio,  i.e.  the  ratio 
°f  lifetime  over  the  transit  time.  Let  us  say  the  ratio  is 
N.  This  phenomenon  results  in  a decrease  of  the  spectrum 
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by  a factor  N2 . When  we  express  the  spectrum  in  terms  of 
equivalent  noise  current,  an  1/N-dependence  results. 

Can  we  explain  the  measured  noise  spectrum  based  on 
this  theory?  If  we  look  at  the  noise  spectrum  between  100 
Hz  and  200  MHz,  we  observe  1/f  noise  like  behavior  with 
several  small  bumps.  Since  the  trapping/detrapping  process 
associated  with  the  quantum-wells  is  quite  fast  due  to  the 
large  capture  cross-sections,  cutoff-frequencies  lie  in 
the  very  high  frequency  domain.  How  do  we  relate  to  this 
low-frequency  Lorentzian  spectrum  then?  We  assume  that 
slow  trapping/detrapping  processes  exist  via  the  traps 
located  in  the  heterostructure  interfaces.  To  have  a 1/f 
noise  like  spectrum  extending  up  to  100  MHz,  we  assume 
with  Tehrani  [71]  that  the  spectrum  is  composed  of 
multiple  Lorentzians  with  parameters  that  are  separated 
more  than  one  decade  from  each  other. 

Thus  we  have  basically  two  kinds  of  traps,  quantum- 
wells  acting  as  giant  traps  and  interface  traps.  The 
latter  is  the  conventional  trap  which  is  caused  by  an 
imperfect  lattice  match  or  by  an  interface  defect  in  the 
heterostructure  junction. 

5.2.  Trapping  and  Detrapping  in  Multiquantum-wells 

Fig.  5.1(a)  shows  the  bandprofile  of  the  MQW 
structures  where  carrier  trapping  and  detrapping  occurs. 
The  thickness  of  the  quantum-well  and  -barrier  are  I2  and 
ll,  respectively.  To  describe  the  dynamics  of  the  hole  and 
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Figure  5.1  (a)  A band  diagram  of  the  multiquantum-well 

heterostructure  layers 
(b)  A simplified  band  diagram  to  model  the 
trapping  and  detrapping  in  quantum  wells. 
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electron  transitions  in  this  system,  a four  energy-level 
system  is  needed.  Finding  an  analytic  solution  for  this 
system  however,  results  in  many  difficulties  due  to  the 
ambipolar  nature  of  the  system  when  we  take  into  account 
both  holes  and  electrons.  To  circumvent  this  situation,  an 
approximation  is  made,  that  is,  when  electron  injection  is 
dominant  as  is  the  case  in  the  devices  under 
investigation,  we  neglect  the  hole  density.  Consequently, 
a simplified  system  consisting  of  two  energy  levels,  i.e., 
the  extended  state  and  the  localized  state  of  the 
conduction  band  shown  in  Fig.  5.1(b),  can  be  used. 

The  basic  equations  which  describe  the  carrier 
transitions  and  flow  are  the  kinetic  equations,  Poisson's 
equation  and  the  current  equation.  For  a two  energy  level 
system,  the  kinetic  equations  read 


In  the  above  equations,  n and  n are  the  electron  density 
in  the  extended  conduction  band  and  the  electron  density 
in  the  localized  conduction  band,  respectively.  The 
Langevin  term,  y represents  the  fluctuations  in  the 
transition  rates  between  the  extended  conduction  band  and 


dn 

dt  “ gn  - rn  + 


(5.1) 


(5.2) 
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the  localized  conduction  band.  The  current  equation  for 
the  electrons,  neglecting  diffusion,  is 


If  we  include  the  displacement  current,  the  total  device 
current  J becomes  equal  to 


where  we  neglect  the  hole  current  as  we  mentioned  earlier. 

To  simplify  the  problem,  we  assume  a quantum-well  to 
behave  as  a giant  acceptor-type  trap.  Thus,  the  quantum- 
well  is  neutral  when  it  is  empty,  whereas  it  is 


the  one-dimensional  Poisson  equation  is  expressed  as 


equal  to  fNA,  and  f is  a trap  filling  factor,  density  of 
the  acceptor-like  traps.  The  concept  of  trap  density  in 
this  context  will  be  discussed  in  section  4 of  this 
chapter . 

For  the  generation  and  recombination  terms  we  have 


Jn  = qjinnE 


(5.3) 


(5.4) 


negatively  charged  when  it  is  filled.  In  this  system, 


(5.5) 


where  E represents  the  electric  field.  The  density  n is 


A 


A 


(5.6) 

(5.7) 


Pn  — ©nn  — Cnnin 


rn  = Cnn<NA  " n) 
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where  Cn  is  a capture  coefficient,  ni  is  the  Shockley 
parameter,  equal  to  the  density  of  electrons  in  the 
conduction  band,  when  the  Fermi  level  coincides  with  the 
trap  level. 

Since  the  total  current  is  solenoidal  Eq.  (5.1)  is 
redundant  and  can  be  dismissed. 

By  substituting  equations  (5.6)  and  (5.7)  into  Eq. 
(5.2)  and  linearization,  we  obtain  the  following  A.C. 
equation  for  the  trapped  electron  density  in  the  quantum- 
wells 


dAn 

dt 


= -A(gn  - rn)  - Y 

= An  (Na  - n0)Cn  - An(ni  + n0)Cn 


Y 


(5.8) 


Now  we  define  several  time  constants 


-1 

Xln 

— Cn  (N^  nQ) 

(5.9) 

-1 

*2n 

= Cn(ni  + nQ) 

(5.10) 

-1 

-l  -l 

(5.11) 

= %n  + T2n 

With  three  definitions,  wee  rewrite  equation  (5.8) 


dAn  _ An  An 

^ ^ln  T2n  ^ 


(5.12) 
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Summarizing,  we  have  three  equations  describing  our 
system.  If  we  express  these  equations  in  a simplified  form 


Ci  = f(J,n,E) 
C2  = f (n,n,E) 
C3  = f(n,n,y) 


(5.4a) 

(5.5a) 

(5.12a) 


where  Ci,  C2,  and  C3  are  constants  and  f denotes  a 
function.  We  have  three  equations  with  five  variables.  Our 
goal  is  to  reduce  these  three  equations  to  one  equation 
which  has  only  three  variables,  i.e.  J,  E,  and  y.  Assuming 
A.C.  short  conditions  (J*Edx  = 0)  or  A.C.  open  circuited 

conditions  (J  = 0),  where  E and  J denote  the  A.C.  Fourier 
components  of  E(t)  and  J(t),  respectively.  In  this  way,  we 
can  obtain  the  current  or  voltage  fluctuations  in  terms  of 
y.  Now  we  return  to  the  reduction  task. 

Fourier  transform  of  Eq.  (5.12)  yields 


(ico  + 


7 


(5.13) 


Upon  substituting  Eq.  (5.12)  into  the  Fourier  transformed 
Poisson  equation,  we  obtain  the  following  first  order 
differential  equation  for  n. 


e_  (l+ioyu2n)  .Jh  . d|  Th  ~ 

q (1+iOJTh)  t2n  dx  (1+iCOTh)  1 


(5.14) 
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To  determine  the  current  fluctuation,  we  return  to  the 
A.C.  current  equation  given  by 

3AE 

Aj  = qp.nQAE  + q|lE0An  + £ (5.15) 

A Fourier  transform  of  this  equation  gives 


J = (qp.n0  + i(D£)E  + qjiE0n 


(5.16) 


Using  the  definition  of  the  dielectric  relaxation  time 


-l  _ q^n0 
G £ 

we  rewrite  Eq.  (5.16) 


(5.17) 


7 = (io)  + 7")E  + -^Vn  (5.18) 

c Iq  c 

where  |iE0  is  replaced  by  vn  which  is  the  electron  drift 
velocity.  By  substituting  Eq.  (5.14)  into  the  above 
equation,  we  obtain  the  following  differential  equation 
for  E . 


dE  _ T2n(l+iCQTh)  (l+i(OTh)  g = _ (l+iC0Th)T2n  g.  + qT2n  ~ 
dx  (l+iOJT2n)%^Vn  ~ £(l+iC0T2n)ThvnJ  + £(l+iOJT2n)Y 

(5.19) 

For  simplicity,  this  equation  can  be  written  in  the 
following  form. 
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(5.19a) 


where 


T2n(l+iQXh) 

(l+i(OT2n)ThTi2Vn 


(5.20) 


(l+i(OTh)^2n  -x.  + ffl2n  ~ 
8(l+iC0T2n)^Vn  e(l+i(OT2n)y 


(5.21) 


The  solution  of  this  differential  equation  which  satisfies 
the  boundary  condition  E = 0 at  x = 0 is  given  by 


Since  the  values  of  P(x" ) and  Q(x')  vary  periodically  in 
the  multiquantum-well  structures,  we  can  write  the  above 
integral  as  a sum  of  sub-integrals.  In  other  words,  we 
need  to  evaluate  a set  of  discrete  integrals,  instead  of 
one  continuous  integral.  To  determine  the  electric  field 
at  each  interface,  we  apply  the  boundary  conditions  shown 
in  Fig.  5 . 1 (b) . 


X 


E (x) 


(5.22) 


J 

o 


E (0)  = 0 
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XI 


Xl 


E (xi)  = e o 


- Jpdx" 


Jpdx" 

Q (x ' ) e o dx ' 


J 

o 


= p-Pll 


i I 


01 


where  Iqj_  represents  the  integral.  Similarly,  at  x = x2. 


E(X2) 


e 


x' 

Jpdx" 

Q (x' ) eo  dx ' 


J 

o 


= p-(Pili+P2l2) 


Uoi  + 


12- 


Note  that  the  main  integral  is  split  into  two  sub- 
integrals, ioi  and  1 12 , because  the  integration  over  x' 
extends  over  two  sections.  In  a similar  fashion,  we  find 
the  electric  field  at  each  interface  X3,  X4,  etc 


E(x3) 


e 


x3 


x’ 

Jpdx" 

Q(x' )e0  dx' 


J 

o 


= e“ (2pi1i+p212) 


<*oi 


+ I12  + I23^ 
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Jpdx" 

Q(x' )e0  dx' 


x 


e-(2P1l1+2P2l2) 


^01  + I12  + x23  + I34^ 


Sub-integrals  101/  Ii2/  123/  etc  are  evaluated  in  the 
Appendix  B.  We  rewrite  the  above  given  expressions  for 
electric  field  by  substituting  the  expressions  of  Ioi/ 
112/  123/  etc. 

E ( 0 ) = 0 


E(x3)  = ^[-e-<2pl1l+p2l2)4e-(pl1l+E‘2l2))^-Pill+1] 

+ ||[  _e-(pl1l+p2l2)4e-Pili] 

E(x4)  = ^[-e-(2Plll+2P2l2)+e-(Pil1+2P2l2)^-(Pili+P2l2)4e-P2l2] 
+ a[^-(P1ll+2P2l2)^-(Pili+P2l2)^-P2l2  +1, 


= 2L[-e-(pill+p2l2)4e-P2l2,+  £-[_e-<Pili+P2l2)+iJ 


Ill 


E(x  ) = Ql[_e-(3Pili+2P2l2)^-(2P1l1+2P2l2)_e-(2P1li+P2l2) 

5 Pi 

4e"^pi1i+P2l2)-e_pi1i  +1] 

+ ^-[-e“  (2Pili+2P2l2)  4e"  (2P111+P212)  _g-  (P111+P212)  ^^1^ 


n,  ,N+1_  , N-l_  , . .N— 1 , .N-l^  , . 

?ir-p_  + e-(— Pili^-^2l2) 


E(xN)  = ~[-e  v 2 


Pi' 


- e-<I!rPi:Ll^P212)-te-(Plll+E'2l2)-e-Plll  +1] 


Oo  /N— 1„  , ,N-1„  , . ,N-1„  , ,N-2^  . 

f— p ( o Pi 3-1  o P2I2)  4p  ( 9 P]_3]d  9 P2I2 ) 


P2 


-e“  (^Pili+^212 ) -e“  (P1I1+P2I2 ) le"15!1!  ] 


if  N is  an  odd  number.  If  N is  an  even  number,  we  obtain 


n,  ,N+1_  , N-l_  , , .N-L  , N-l  , , 
= ?ir-P~(~pi1i4'T‘f2l2)+  e-(— PiliP^-P2l2) 


E(xn)  = — [~e  v 2 


Pi1 


- e-(^P1ll^P2l2)^-(P1l1+P2l2)4e-P2l2  ] 


ry,  ,N„  n ,N_  , . .N-2^  n ,N_  . . 

fp^[-e  ?1^"1+2P2^"2  ) +e  ^ 2 Pl-*-l+2P2l2) 


-0""  (~2~P111"^2~P212)  +e-(Pili+P2l2)_^-P2l2+  1] 


As  shown  in  Fig.  5.1b,  the  last  section  of  the 
multiquantum-well  layers  has  a small  band-gap.  This 
implies  that  N is  an  odd  number  in  our  situation. 


112 


Note  that  the  rightmost  term  in  each  bracket  that 
has  a coefficient  Q1/P1  is  equal  to  unity  when  N is  an  odd 
number.  If  N is  an  odd  number,  the  rightmost  term  in  each 
bracket  which  has  a coefficient  Q2/P2  is  equal  to  unity. 
Hence  summation  of  the  electric  field  at  each  interface 
results  in  the  following  expression 


2> 


(xi)—  1+0+1+0+ 

+ ^j[0+l+0+l+ 


+0+l+exponential  terms] 
+l+0+exponential  terms] 


(5.23) 

where  the  contribution  of  the  exponential  terms  is 
negligible  with  respect  to  unity  because  of  the  following 
reasons.  Firstly,  each  exponential  term  has  an  alternating 
sign,  thus  reducing  the  contribution  to  the  overall 
summation.  Secondly,  each  exponential  term  becomes 
smaller  as  N increases.  Therefore  we  can  rewrite  Eq. 
(5.23) 

N 

Y~  N Qi  Q? 

E(xi)=  (5.24) 

The  voltage  across  the  depletion  region  is  given  by  the 
line  integral  of  the  electric  field,  i.e.. 
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L 


0 


For  the  multiquantum-well  structures,  we  get 


(5.25) 


V 


N 

{ SE(xi)Axi 

i=odd  num. 


+ 


N 

5^E(xi)Axi  } 

i=even  num. 


Nil  Ql  Nl2  Q2 

2 Pi  2 P2 

To  evaluate  Eq. (5.26)  we  use  Eq. (5.20) 
modified  Eq.  (5.21)  as  follows 


(5.26) 

for  Pi,  P2  and 


(l+iC0Th)X2n  ~ 
£ (1+i  (0l2n ) Th  vn 


(5.27) 


where  we  put  y equal  to  zero,  since  we  are  interested  in 
finding  the  small-signal  impedance.  Using  Eq.  (5.26)- 
(5.27).  The  impedance  is  found  from 


7 _ v_  _ NiiQi  ni2  Q2 

JA  = 2 pl  2 p2 

Let  Li=Nli/2  and  L2=Nl2/2,  then 


(5.28) 


LlQl  L2Q^ 
PlA  + P2A 


(6.29) 


Upon  substituting  Eqs.  (5.20)  and  (5.27)  into  (5.29) 


Z = 


Li 


£lA(iC0P— ) 
tia 


+ 


L2 


1 

e2A(lO)4— ) 
T2£2 
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Ll  + L2 

aiA(l+icoxm)  a2A(l+icoT2n) 


(5.30) 


where  we  used  the  definition  of  the  dielectric  relaxation 
time  i.e.,  Tq  = e/a.  To  illustrate  the  physical  meaning  of 
Eq.  (5.30),  we  use  the  following  definitions 


L , 8A 

R = — and  C = -r- 

GA  L 

Consequently,  Eq. (5.30)  yields 
R1  R2 

Z — -}-  — 

l + iO)RiCi  l + i(0R2C2 


(5.31) 


As  Eq. (5.31)  indicates,  the  overall  impedance  of  the 
multiquantum-well  structures  is  equal  to  the  series 
connection  of  two  kinds  of  impedances,  while  the  one  with 
subscript  1 corresponds  to  the  small  band-gap  section, 
while  the  other  with  subscript  2 is  the  impedance  of  the 
large  band-gap  section.  Each  impedance  can  be  represented 
tas  a resistor  and  a capacitor  in  parallel. 

For  later  use,  we  define  an  impedance  correction 
factor,  fz  as  follows 


Li  e2Tm(l+icoT2fl) 

L2  eiT2Q(l  + itt>T2fi) 


(5.32) 


Then  Eq.  (5.30)  can  be  rewritten  as 

2 _ f ^20^2 

2 e2A(l  + iCOT2Q) 


(5.33) 


In  our  case  Li  = L2  = L/2.  Then,  Eq.  (5.33)  yields 


Z 


l2Ql 

z 2e2A(l+icoT2n) 


(5.34) 


Eq.  (5.34)  will  be  used  later  to  determine  the  current 
fluctuation  noise. 

The  transport  equation  (5.19)  that  we  have  used  so 
far  can  yield  the  noise  also.  Previously,  we  neglected  the 
fluctuation  terms  to  find  the  impedance.  Now  we  put  J = 0 
in  Eq.  (5.21)  to  determine  the  terminal  noise  in  terms  of 
open-circuited  voltage  fluctuations.  We  rewrite  Eq.  (5.21) 
with  J = 0 


Q = 


~ 

e(l+ian&n) 


(5.21a) 


For  the  open-circuit  voltage  we  find 

L 

V (L,  (0)  = - E(x,CG)dx  (5.35) 

<? 

where  the  electric  field  expression  contains  the 
fluctuations  implicitly.  The  spectral  density  of  the  open- 
circuit  voltage  fluctuation  can  be  found  from 
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where  < > denotes  average  value  and  T is  a given  time 
interval.  To  evaluate  Eq.  (5.33),  we  split  the  integral 
into  N sections  as  we  did  for  the  impedance.  Eq. (5.36)  can 
be  rewritten  as 


SV  = 7f<  t 


[/s’ 


Xl  x2 

J*E(x,  (0)dx+  j*E(x, 
X2 


x3 


XN 


CO)dx+ Je(x,  C0)dx+- • • + J*E(x,  co)dx] 

Xfj-i 


^2 

X3 

Ji 


[ |E*(x',co)dx'  + E*(x',co)dx'  + E (x',co)dx'  + - • •+  I E (x ' , co)dx'  ] > 

Xl  X2  *N~1 


Xn 


Xl 

r xi 


T [ 

J 

0 

x3 

(*  x3 


J*£(x, 
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J 
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J^E(x, 


co)  E*  (x' , co) dxdx' + 


J 

0 


x2 


(5.37) 

The  evaluation  of  the  double  integrals  in  the  above 
equation  is  given  in  the  Appendix  C.  Note  that  the  open- 
circuit  expression  for  the  electric  field  E(x,co)  has 
fluctuation  terms  only.  Hence  the  product  of  E(x,co)  with 
E * ( x , co)  yields  the  product  of  fluctuation  terms  i.e., 
Y (x,  co)y*  (x ' , co)  . The  average  value  of  this  product  is 
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proportional  to  a spectral  density  and  is  given  by  g-r 
noise  theory.  If  we  write  this 


Since  Sy  contains  delta  function,  the  cross  terms  in  the 
last  line  of  Eq.(5.37)  become  zero.  This  means  that  g-r 
noise  sources  in  each  quantum-well  are  statistically 
independent  of  each  other.  Also  note  that  the  double 
integrals  become  single  integrals  because  of  the  delta 
function  in  the  integrand. 

Using  a result  we  obtained  from  Appendix  C 
xn 


(5.38) 


2 

T 


Ji 


E(x, (0)E*(x’,(0)dxdx'~  0;  If  N is  an  odd  number 


We  rewrite  Eq. (5.37)  and  find 


E (x,  co ) E*  (x ' , (0 ) dx  dx 
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Nl?  .Q2Q2* 
T P2P2 


(5.39) 


For  the  structures  where  li  = I2  = 1/  Eq.  (5.39)  can  be 
rewritten  as 


Svh 


Nl  Q2Q2* 
T P2P2 


To  evaluate  Eq.  (5.39a),  we  use  Eq.  (5.20), 
(5.38) , and  find 


(5.39a) 
(5.21)  and 


Svn  = Nl  — 

82  (l+C02Xn)  (l+CO2^) 


(5.40) 


Concerning  the  fluctuation  terms  in  Eq. (5.46),  we  assume 
that  the  random  emptying  and  filling  of  traps  results  in 
full  shot  noise,  allowing  us  to  write  the  noise  spectrum 
as  follows 


4nQ 

Alin 


(5.41) 


Using  this  expression,  Eq.  (5.40)  can  be  rewritten  as 


Svn  = Nl  ~ 

£2  (l+CO2^) 

4Nl|q2n0vnin 


A82  (1+(02X2) 


,4n0 

AX 


) 


(5.42) 
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where  Xfn  is  the  flight  time  for  the  electrons  between  the 
adjacent  quantum-wells.  In  Eq.  (5.42),  we  used  vn  = l/ifn. 
The  spectral  density  of  the  current  fluctuations  can  be 
found  with  the  help  of  Eq.  (5.41) 


4Nl2q2n0vnXQ  4e2A2  ( 1+C02Xq) 

_ ^vn  _ X]_n  1 

Sln  “ |Z|2  “ A£2(1+(o2t2)  <~)  $2i2t2  ’ if 

= (5.43) 

Note  that  in  Eq.  (5.32)  the  impedance  correction  factor  fz 
is  approximately  2 at  low  frequencies.  Using  the  value  of 
fz  and  I = qn0vnA,  Eq.  (5.43)  can  be  rewritten  as 


Si 


n 


iSl/IlUx 

N ' Xf 


(5.44) 


From  this  result,  we  note  that  the  current  noise  level 
critically  depends  on  the  number  N of  heterostructure 
layers  and  on  the  ratio  (Xin/Tf)  . The  bias  dependence  of 
the  noise  level  can  be  explained  in  terms  of  these  two 
parameters.  If  the  quantum-wells  are  depleted  by  applying 
the  high  bias  voltage,  then  xin  = NXf  and  full-shot  noise 
is  obtained.  At  low  bias,  however,  the  time  constant  ratio 
becomes  one  in  the  undepleted  quantum-well  layers  and  a 
strong  suppression  of  the  noise  will  occur  due  to  a factor 
of  1/N,  where  N is  the  number  of  undepleted  quantum-wells . 
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As  we  will  see  in  chapter  VI,  as  reverse  bias  increases, 
the  quantum-well  layers  in  the  i-region  are  depleted 
gradually  i.e.,  one  quantum-well  at  a time.  Thus, 
reduction  in  the  number  N of  the  undepleted  quantum-wells 
results  in  the  increase  in  the  noise  level. 

5.3.  Trapping  in  Heteroiunction  Interface  States 

As  mentioned  in  the  introduction,  the  development  of 
the  theory  presented  in  this  section  was  prompted  by  the 
frequency  dependence  of  the  noise  spectrum  in  the  mid- 
frequency range.  We  assume  the  levels  of  traps  in  the 
forbidden  band-gap  as  shown  in  Fig.  5.2a.  The  difference 
with  the  conventional-type  traps  is  that  the  traps  in  this 
case  are  located  at  discrete  positions,  i.e.,  only  at  the 
heterostructure  interfaces. 

Fig.  5.2a  shows  the  band  diagram  for  the  interface 
trap  model.  In  order  to  contribute  to  an  electrical 
conduction,  the  electrons  in  each  trap  level  should 
transit  to  the  energy  level  of  the  extended  conduction 
band  (Eec) . This  electron  transition  consists  of  a 
transition  between  a trap  energy  level  (Et)  and  the 
localized  conduction  band  (Eic)  , followed  by  a transition 
from  the  E]_c  state  to  the  Eec  state.  The  latter  is  the 
quantum-well  detrapping  and  detrapping  process  discussed 
in  the  previous  section.  Since  the  quantum-well  trapping 
and  detrapping  is  a very  fast  process,  the  above  two-step 
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Interface  States 


4 

3 

2 

l 


(a) 


(b) 


Fig  5.2  (a)  Hetero junction  interface  traps 
(b)  simplified  interface-trap  model 
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process  can  be  approximated  by  a one-step  process,  namely, 
the  Et-Eic  transition  is  assumed  to  be  identical  to  the 
ET“Eec  transition.  This  transition  is  shown  in  Fig. 
5.2 (b) . 

With  the  carrier  transitions  and  flow  as  indicated  in 
Fig.  5.1(b),  the  kinetic  equations  become 


dn  1 dJn 

dt  = gn  - rn  + + Y 


(5.45) 


dn 

— = -gn  + rn  - y 


(5.46) 


In  the  above  equations,  n and  n are  the  free  carrier 
density  and  the  trapped  carrier  density,  respectively,  y 
represents  the  fluctuations  in  the  transition  rates.  If  we 
neglect  diffusion  current  and  include  the  displacement 
current 


J - 


_ 3e 
qn|inE  + fr^-. 


Poisson's  equation  becomes 


(5.47) 


dE 

dx 


q , a 

j(n  + n - Nd), 


(5.48) 


where  denotes  the  impurity  density. 

For  the  generation  and  recombination  terms  we  have 


gn  - enn  = Cnnin 


(5.49) 


rn  = Cnn(Nd  - n) 


(5.50) 
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where  Cn  is  the  capture  coefficient,  and  ni  is  the  density 
of  electrons  in  the  conduction  band  when  the  Fermi  level 
coincides  with  the  trapping  level . 

By  substituting  equations  (5.49)  - (5.50)  into  Eq. 
(5.46)  and  linearization,  we  have  the  following  A.C. 
equation  for  the  trapped  electron  density  in  the  interface 
states 


= An  (Nt  - n)Cn  - An  (ni  + n0)Cn  - 7 


An  An 
Ti  12  ” Y 

(5.51) 

with 

-1 

T1 

= Cn(Nd  - h) 

(5.52) 

-l 

*2 

= Cn(ni  + nG) 

(5.53) 

Let 

I"1 

-l  -l 

Xl  + T 2 

(5.54) 

We  apply  Fourier  analysis  on  the  linearized  A. C. equations 
of  (5.47),  (5.48)  and  (5.51).  With  some  algebra,  we  obtain 

the  following  first  order  differential  equation 

dE  X2  (1+iftre)  (l+ia)Tfl)  ^ (l+ia>x)T2  ~ gi2  ~ 
dx  (l+i(0T2)Txgvn  E e(l+icar2)tvnJ  + e(l+icai2)^ 

We  rewrite  Eq. (5,55)  in  the  following  form 


dE 

dx 


+ PE 


Q 


(5.55a) 
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where 

= T2  (1+i  cor)  (1+icoTh) 

( 1+i  0JU2 ) 


(5.56) 


(1+i  cot)  t2  ~ gr2  ~ 

e(l+icoT2)xvnJ  e(l+icoT2)  ’’ 


(5.57) 


Solving  Eq.  (5.55),  we  obtain  for  the  spectral  density  of 
the  voltage  fluctuations 

q _ N1  Q2Q2* 

Sv  = T l p2p*  > Ky  (5.58) 

where  subscript  ' 2 ' denotes  the  value  of  the  even-numbered 
section.  Substitution  of  Eqs.(5.56)  and  (5.57)  and  (5.58) 
yields 

4 N 1 q2  n0v^T2 

Sy  = (5.59) 

Ae2(l+co2x^)  (l+CG2T2)Tl 

The  current  spectral  density  is  found  with  the  help  of  the 
impedance 

4Nlq2n0v^T2TQ  e2A2  (1+(02^) 

Si  - — = . 

*Z*2  AS2  (1+C02X^)  (1+C02T2)Ti  N212T^ 


4qI  tLXl  (1+C02T2) 


(5.60) 


Eq.  (5.60)  is  a well  known  g-r  noise  spectrum  result  and 
shows  the  Lorentzian-type  frequency  response. 
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So  far  we  have  considered  only  a single  trap  level . 
For  the  multilevel  trap  system,  we  follow  the  same 
procedure  as  Tehrani  [71]  did.  That  is,  we  assume  a linear 
superposition  of  multiple  Lorentzian  components,  each  with 
different  trap  parameters.  Thus  we  get  the  following 
spectrum 


Si  = 4ql  2 — rs  ~ 

i tL  Tii  ( 1+C02T2) 


(5.61) 


where 


- Cni  (Nti  nj_) 


-l 

^2i  - Cn  (nil  + no) 


(5.62) 

(5.63) 

(5.64) 


5.4.  Computer  Calculations  of  the  Trapping  Noise 
It  is  difficult  to  present  a quantitative  analysis  of 
the  trapping  noise  of  the  multiquantum-well  structures 
because  of  the  many  unknown  parameters.  For  this  reason  we 
can  only  hope  to  carry  out  qualitative  calculations  based 
on  the  limited  information  available  from  the  experimental 
results . 

We  begin  with  a calculation  of  the  capture  cross 
section  of  the  quantum-well  traps.  Unlike  the  capture 
cross  section  of  an  ordinary  trap  center  (see  Fig.  5.3) 
which  has  a typical  value  of  10-20  m2,  a heterostructure 
quantum-well  has  a very  large  cross  section,  i.e.  ~ 10-8 
m2,  because  the  capture  cross  section  is  associated  with 
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the  total  area  of  the  heterostructure  interface.  If  we  use 

this  cross  section  value,  we  obtain  a very  small, 

unrealistic  value  for  our  time  constant.  Therefore,  we 

need  to  modify  the  definition  of  capture  cross  section  in 

the  context  of  quantum-well  structures. 

For  free  carriers  to  be  captured  by  a trap,  they  need 

to  lose  energy  through  a scattering  process.  According  to 

Whiteside  [72],  the  scattering  mean  free  path  is  of  the 
o 

order  of  10  A for  energies  of  about  1.0  ev.  Then  we  are 

o 

allowed  to  assume  that  at  every  10  A,  free  carriers 
experience  scattering,  loose  energy  and  can  be  captured  by 
the  traps.  This  idea  introduces  the  following 
approximation.  Instead  of  modifying  the  capture  cross 

° o 

section,  we  use  the  standard  value  of  (10  A) ^ and  we 

o 

assume  that  the  traps  are  spaced  10  A apart  forming  a 
lattice  as  shown  in  Fig.  5.4.  Then  we  can  calculate  an 
effective  trap  density 

N = JUgfcer  of  traps  = 1 _ 1027m-3  (5.65) 

unit  volume 

(10A)  j 

This  value  leads  to  values  of  Xin  and  X2n  of 


tin  — 


CnNA 


= (5cvthNA)  1 ~ 3.8  x 1013 


(5.66) 


As  shown  in  Eq.  (5.43)  the  evaluation  of  the  noise  spectra 
depends  critically  on  the  time  constant  Xin.  Calculation 
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Figure  5.3  The  apture  cross  section 
of  ordinary  trap  centers. 


Figure  5.4  The  effectice  capture  cross  section 
of  heterostructure  quantum-wells. 
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of  the  current  noise  spectrum  associated  with  quantum-well 
trapping  can  be  performed  by  using  Eq.(5.31),  (5.43)  and 
the  parameters  presented  in  Table  5.1. 


Table  5.1 

Parameters  used  for  the  calculation 
of  the  quantum-well  trapping  noise 


I,  device  current 

1 HA 

L,  width  of  the  NQW  region 

2-10-6  m 

N,  number  of  heterostructure  layers 

200  (for  UM505) 

A,  cross  section  area  of  the  device 

3.14  -10-8  m2 

vn,  electron  velocity 

1-105  m/sec 

Cnw,  capture  rate  of  the  quantum-well  traps 

2.6-10-15  m/sec 

Nt,  density  of  quantum-well  traps  in  Ec 

1-10 27  m-3 

tL,  transit  time  across  the  MCW  region 

2 -10_11sec 

if,  flight  time  between  each  quantum-veil 

2- io -13 sec 

lb,  dielectric  relaxation  time 

2 -10_6sec 

£,  relative  dielectric  constant 

12 

£or  free  space  dielectric  constant 

8. 85- 10  “12  f/m 

The  result  of  this  calculation  indicates  that  the 
spectrum  is  flat  up  to  tens  of  GHz.  As  expected  in  the 
introduction,  the  magnitude  of  the  low-frequency  spectrum 
plateau  value  is  reduced  by  a factor  of  1/N,  where  N is 
the  number  of  the  heterostructure  layers  in  the 
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multiquantum-well  device.  In  other  words,  this  distinctive 
sub-shot  noise  phenomenon  is  strong  evidence  for  the 
occurrence  of  the  trapping  and  detrapping  in  the  quantum- 
wells  . 

Next  we  focus  on  Eq.  (5.61).  To  explain  the  measured 
noise  spectra  which  show  approximately  1/f  noise  like 
behavior,  a superposition  of  several  Lorentzians  with 
different  time  constants  was  considered  by  Tehrani  [71] . 
As  indicated  by  van  Rheenen  [73]  superposition  to  1/f  is 
possible  only  when  the  time  constants  are  a factor  of  ten 
or  more  apart.  Then,  we  are  allowed  to  treat  a trap  level 
in  the  forbidden  gap  separately  and  single  trap  level 
theory  applies  to  the  multi-level  trap  situation. 

Eq. (5.61)  shows  that  each  Lorentzian  is  characterized 
by  a plateau  value  and  a time  constant  whose  reciprocal 
value  is  the  cutoff  frequency.  To  have  a 1/f  noise  like 
spectrum,  we  noted  that  plateau  value  of  each  Lorentzian 
should  vary  inverse  proportionally  to  the  frequency.  Since 
the  transit  time  tL  has  a fixed  value,  having  a 
statistical  factor  x/ii  independent  of  the  degree  of  trap 
occupancy  is  able  to  realize  latter  condition.  The  time 
constants  T and  X\  are  known  to  be  a function  of  the 
following  trap  parameters:  activation  energy,  carrier  and 
trap  densities,  spin  degeneracy,  and  capture  cross 
section.  To  evaluate  Eq. (5.62)  and  (5.63),  we  present  the 
following  expressions. 
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The  average  value  of  the  trapped  carrier  density  is 
given  by 


n0 


Nt 


1 Et-Ef 

lHg®xp  'n*-* 


(5.67) 


where  g denotes  spin  degeneracy.  The  average  free  carrier 
density  and  the  Shockely-Read  parameter  is  expressed  by 


nQ  = Nc  exp  ( 


EF~EC, 
kt  ’ 


(5.68) 


ni 


Et~E 

kt 


(5.69) 


respectively. 

To  explain  the  measured  noise  spectrum,  we  need  to 
extract  an  approximate  value  for  the  statistical  factor 
x/xi  from  Fig.  5.5.  If  we  substitute  the  value  of  the 
spectral  density  at  the  corner  frequency  into  Eq.(5.60), 
we  obtain  the  following  condition 


— = 10-3  ~ 10-5  (5.70) 

Tl 


Using  Eq.  (5.62)  and  (5.63),  we  find 


12  ^ _L  = Nt-n0 

xi  xi  ni+nQ 

If  nQ  is  negligibly  small  compared  to  Nt  and 


(5.71) 
if  nQ  is 


smaller  than  ni,  the  statistical  factor  is  about  constant. 
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Then,  by  comparing  Eq.(5.68)  and  (5.69),  we  note  that  the 
Fermi  energy  level  must  be  located  below  the  trap  energy 
level.  Since  our  multiquantum-well  structures  are  undoped, 
we  assume  that  the  Fermi  energy  level  lies  around  mid-gap. 
Next  we  assume  the  presence  of  several  discrete  trap 
levels  located  between  the  Fermi  level  and  the  conduction 
band.  We  label  the  trap  levels  as  shown  in  Table  5.2. 
Using  these  parameters,  we  obtained  a 1/f  noise  like 
spectrum.  If  we  combine  this  result  and  the  result  of  the 
high-frequency  spectrum,  we  are  able  to  explain  the  entire 
spectrum  ranging  from  1 KHz  up  to  100  GHz  (see  Fig.  5.5). 


Table  5.2 

Trap  parameters  used  for  the 
calculation  of  the  interface  trapping  noise 


Trap  # 

Nt  (m-3) 

Eact (ev) 

Cn (m3/sec) 

x (sec) 

fc(Hz) 

1 

1014 

0.318 

10-15 

10“3 

103 

2 

1015 

0.259 

10-15 

10"4 

104 

3 

1016 

0.199 

10-15 

10"5 

105 

4 

IOI7 

0.139 

10-15 

10"6 

106 

5 

lOl8 

0.079 

10-15 

lO-7 

107 

6 

1019 

0.019 

10-15 

10"8 

108 
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Figure  5.5  Noise  spectrum  obtained  by 
the  computer  simulation 


CHAPTER  VI 

EXPERIMENTS  ON  MULT I QUANTUM-WELL  (MQW) 
HETEROSTRUCTURE  P-I-N  PHOTODIODES 

6.1.  Introduction 

An  experimental  study  on  heterostructure  p-i-n 
photodiodes  is  discussed  in  this  chapter.  The  emphasis  is 
placed  on  the  noise  measurements  and  their  analysis. 
Together  with  noise  measurements,  conventional 
measurements  such  as  I-V,  C-V,  impedance  (S-parameters) 
should  be  performed,  since  these  measurements  provide 
information  about  the  steady-state  operating  condition  of 
the  device.  Researchers  have  been  using  noise  measurement 
as  a powerful  tool  in  studying  the  microscopic  behavior  of 
an  electronic  device  and  its  physics. 

Our  initial  objective  of  this  study  was  to  measure 
the  impact-ionization  coefficient  ratios  on  the 
heterostructure  MQW-APDs  using  the  noise  measurement 
method.  When  it  turned  out  that  we  could  not  obtain  a 
high  enough  electric  field  to  initiate  the  impact 
ionizations  in  the  InGaAs  based  structures  due  to  the  high 
tunneling  current  at  room  temperature,  we  discovered  that 
some  multiquantum-well  structures  exhibited  distinctive 
sub-shot  noise  and  that  this  sub-shot  noise  is  dependent 
on  both  the  structure  and  the  applied  bias.  So  it  was 
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decided  that  sub-shot  noise  needed  to  be  studied  instead. 
Measurements  on  three  different  quantum-well  structures 
were  performed  and  their  findings  will  be  presented  in 
this  chapter. 

In  the  first  section,  we  will  discuss  the  device 
structures.  In  subsequent  sections,  the  experimental  set- 
up, the  experimental  procedures  and  the  results  will  be 
presented  with  a discussion. 

6.2.  Device  Description 

A heterostructure  APD  in  its  simplest  form  is  a p+-i- 
n+  diode.  Instead  of  an  i-region,  an  undoped  MQW  layer  is 
sandwiched  between  p+  and  n+  layers.  The  configuration  of 
the  device  used  in  the  experiments  is  shown  in  Fig.  6.1. 
The  multiquantum-well  layer  and  absorption  layers  were 
grown  on  an  n+  InP  substrate  by  molecular  beam  epitaxy 
(MBE)  technique  at  the  University  of  Michigan. 

For  device  UM505,  a 2 (im-thick  n+  Ino.53Gao.47As  layer 
for  photo-absorption,  a 1.8  Jim-thick  undoped 
In0.53Ga0 . 47As/In0 . 52AI0 . 48As  MQW  region  and  another  2 fim- 
thick  p+_Ino . 53Gao . 47AS  layer  for  top  illumination  were 
subsequently  grown.  Three  different  devices  with  different 
MQW  periods  and  dimensions  were  used  in  the  experiments. 
Detailed  information  of  these  devices  is  given  in  Table 
6.1. 
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Bonding  Wire 


■ 


* 


Optical  Fiber 


Figure  6.1  Configuration  of  a multiquantum-well 
heterostructure  p-i-n  photodiode 
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Table  6.1 


Device  No.  Composition 

No.  of 
MQW  periods 

Barrier 

Width 

Well 

Width 

UM504  Ino.53Gdo.47As/  lno.52Alo.4gAs 

180 

30  & 

30  \ 

UM505  Ino.53Gao.47As/Ino.52Alo.48As 

100 

90  & 

90  \ 

UM516  Ino.53Gao.47As/Ino.52Alo.48As 

25 

500  & 

500  A 

A window  is  provided  in 

the  metal 

film  on 

the  top 

p+-InGaAs  layer  for  illumination.  This  top  illumination 
gives  rise  to  a photo-generation  of  electrons  and  holes. 
The  holes  are  collected  by  the  top  electrode  and  excess 
electrons  will  be  injected  into  the  multiquantum-well 
layers . 

When  this  photodiode  is  operated  at  high  reverse 
bias,  most  of  the  voltage  drop  is  across  the  undoped  MQW 
region,  resulting  in  high  electric  fields,  and  thus 
providing  the  possibility  for  multiplication.  The 
multiplication  process  can  not  take  place  in  the  n+  and 
p+  InGaAs  layers  because  the  voltage  drop  across  the 
heavily  doped  layers  is  negligibly  small.  These  layers, 
however,  consist  of  low  band-gap  material  (Eg  =0.75  eV) 
which  is  suitable  for  photo-absorption  of  infrared  (1.3  mm 
and  1.5  mm)  wavelengths  where  material  dispersion  and 
optical  attenuation  are  minimum  in  silica  based  optical 
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fibers.  The  multiquantum-well  layer  consists  of  two 
different  band-gap  materials,  i.e.,  Ino.53Gao.47As  (Eg  = 
0.75eV),  and  Ino . 52-Alo . 4 8As  (Eg  = 1.47  eV)  . At  the 
heterostructure  interface  the  difference  in  band-gap 
results  in  band-edge  discontinuities.  Recent  measurements 
have  shown  that  70%  of  the  band-gap  difference  lies  in  the 
conduction  band.  Thus  we  assume  a conduction-band 
discontinuity  of  0.5eV  and  a valence-band  discontinuity  of 
0 . 2eV  [74]  for  our  devices.  Since  carrier  transport  is 
perpendicular  to  the  heterostructure  interface,  electrons 
and  holes  experience  different  potential  discontinuities 
and  consequently  carriers  will  lose  or  acquire  different 
amounts  of  kinetic  energy  at  each  interface. 

6.3.  D.C.  Measurements  and  Results 

Steady  state  characterizations  such  as  I-V,  C-V,  and 
impedance  measurements  are  necessary  to  understand  the 
device  operation  and  need  to  be  done  before  the  noise 
characterizations.  Measurements  were  performed  at  room 
temperature  on  the  three  devices  with  different  MQW 
structures . 

We  begin  with  device  UM504  (30  X - 30  X)  . A 
semiconductor  parameter  analyzer  (HP4145A)  was  used  to 
measure  the  reverse  dark  current-voltage.  The  result  is 
shown  in  Fig. 6.2.  At  low  bias  the  reverse  dark  current  is 
dominated  by  generation-recombination  process  in  the  MQW 


Current  (A) 


138 


Reverse  Bias  Voltage 


Figure  6.2  Current -voltage  (I-V)  characteristics. 
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structure  rather  than  by  diffusion  in  the  InGaAs 
absorption  regions,  because  of  the  low  density  of  minority 
carriers  in  this  highly  doped  material.  The  rapid  increase 
in  the  dark  current  above  15  V is  due  to  Zener  tunneling 
in  the  low  band-gap  layers  (Ino.53Gao.47As) . Photocurrent 
was  generated  by  illuminating  the  top  layer  (p+-InGaAs) 
with  a LED  at  A,  = 1.3  |im.  The  reverse  current-voltage 
characteristic  was  measured  with  continuous  bias.  Between 
0 and  7 V the  increase  in  the  photocurrent  is  due  to 
enhanced  injection  of  excess  electrons  into  the  depletion 
region  as  bias  increases.  Above  7 V,  the  photocurrent 
remains  constant  for  a while  and  increases  again  at  higher 
bias  due  to  the  increasing  dark  current  component.  From 
this  I— V characteristic,  we  see  that  this  device  performs 
similarly  to  an  ordinary  bulk  device.  We  did  not  increase 
the  bias  high  enough  to  observe  the  impact  ionization 
because  the  tunneling  current  is  already  too  high. 

Fig.  6.3  shows  the  result  of  the  C-V  measurements. 
The  rapid  decrease  in  the  junction  capacitance  between  0 
and  0.5  V indicates  that  the  undoped  MQW  layers  are  easily 
depleted  by  a small  bias.  Further  decrease  is  due  to  the 
depletion  at  the  edge  of  the  InGaAs  absorption  layers. 
Above  12  V,  the  junction  capacitance  Cj  reaches  a constant 
value  of  2.77  pF  which  is  quite  close  to  the  geometrical 
value  of  2.71  pF. 
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Figure  6.3 
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characteristic 
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The  impedance  was  measured  with  an  S-parameter 
analyzer  (HP  8510A) . The  Smith  chart  plot  of  the  result  is 
shown  in  Fig.  6.4.  To  extract  the  series  resistance  value 
we  used  simulation  software  (Touchstone)  and  obtained  the 
value  of  30-Q.  This  value  together  with  the  junction 
capacitance  will  be  used  in  the  noise  data  analysis. 

Measurements  of  the  I-V  characteristics  of  device 
UM505  (90  X - 90  X)  are  shown  in  Fig.  6.5.  As  in  device 
UM504,  below  15  V,  the  dark  current  is  dominated  by 
generation-recombination  processes  and  diffusion.  With 
increasing  bias,  the  depletion  of  more  quantum-well  layers 
results  in  a rapid  increase  in  the  field  assisted  thermal 
generation  in  the  depleted  layers.  At  extremely  high  bias, 
Zener  tunneling  current  becomes  dominant. 

Note  that  the  photocurrent  is  strongly  dependent  on 
reverse  bias.  Since  the  carrier  generation  by  photo- 
absorption is  independent  of  bias,  this  bias  dependent 
quantum  efficiency  is  attributed  to  the  bias  dependent 
collection  process  of  the  carriers.  In  an  ordinary 
homostructure  p-i-n  photodiode  with  separated  absorption 
layer,  the  carriers  are,  upon  injection  into  the  depletion 
region,  accelerated  by  the  electric  field  and  cross  the 
depletion  region.  Due  to  a lack  of  mobile  carriers  in  this 
region  recombination  is  unimportant.  If  there  is  a 
potential  barrier  between  the  absorption  layer  and 
depletion  region,  however,  as  in  the  heterostructure 
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Figure  6.4  A plot  of  S-parameter  (S,-,  ) 
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Figure  6.5  I-V  characteristics  of  the  device  UM505. 
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photodiode  shown  in  Fig.  6.6,  photocurrent  will  be 
strongly  dependent  on  the  amount  of  carriers  injected  over 
the  first  few  potential  barriers.  Effectively,  these 
barriers  reduce  the  carrier  injection  into  the  depletion 
region.  In  other  words,  an  increase  in  the  bias  will  lead 
to  a reduction  in  the  effective  barrier  height  and 
consequently  to  an  increase  in  the  photocurrent.  In  our 
structure  (UM505) , the  collection  efficiency  increases 
gradually  as  a function  of  bias.  To  understand  this,  we 
present  in  Fig.  6.7  the  successive  band  profiles  with 
increasing  bias.  As  shown  in  this  figure,  an  increase  in 
the  bias  results  in  the  depletion  of  more  quantum-wells 
instead  of  lowering  the  potential  barrier  by  tilting  the 
band-profile.  This  weak  dependence  of  the  barrier  height 
on  the  reverse  bias  accounts  for  the  gradual  increase  in 
the  quantum  efficiency  at  low  bias.  In  the  previous  device 
(UM504)  which  has  a very  thin  potential  barrier,  we  didn't 
observe  this  effect,  because  the  carriers  can  easily 
tunnel  through  the  barriers,  whereas  in  an  InP/InGaAs 
single  hetero junction,  nontransparent  structure,  Forrest 
et  al . [75]  observed  that  unity  collection  efficiency  is 
obtained  rather  abruptly. 

At  higher  bias,  however,  the  potential  barrier 
limiting  the  carrier  injection  is  fully  lowered  and  the 
photocurrent  is  limited  by  carrier  transit  across  the 
undepleted  quantum-well  layers. 
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p+-  InP  inP  n-  InGaAs 


(a) 


Figure  6.6  (a)  Equilibrium  and  (b)  reverse-bias  energy  band 

diagrams  of  an  InGaAs-InP  photodiode  with  separate 
absorption  and  multiplication  regions. 
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(a) 


Figure  6.7  Successive  depletion  of  quantum-well 
layers  from  (a)  low  to  (d)  high  bias 
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C-V  measurements  are  not  available  for  device  UM505  , 
but  the  calculated  value  of  the  geometrical  capacitance 
yields  2.31pF.  From  this  value  and  from  the  result  of  the 
S-parameter  measurement  shown  in  Fig.  6.8,  we 
approximately  obtained  a series  impedance  value  of  10-Q. 

Device  UM516  was  measured  next.  As  shown  in  Fig.  6.9, 
the  I-V  characteristic  demonstrates  the  gradual  increase 
in  the  quantum  efficiency  as  a function  of  reverse  bias. 
The  results  of  the  C-V  measurement,  shown  in  Fig.  6.10, 
also  indicate  the  gradual  depletion  of  the  multiquantum- 
well  layers  over  a wide  bias  range,  i.e.  between  0 and  30 
V.  The  reason  for  the  sudden  increase  in  the  junction 
capacitance  near  zero  bias  is  unknown.  The  measured 
junction  capacitance  Cj  at  high  bias  yields  1.35  pF  which 
is  close  to  the  geometrical  value  of  1.17  pF.  Fig.  6.11 
shows  the  results  of  S-parameter  measurements.  From  this 
data  we  obtained  a series  resistance  value  of  10-Q. 

6.4.  Noise  Measurement  Set-up 
One  of  the  common  difficulties  in  noise  measurements 
is  to  amplify  the  small  noise  signal  to  a measurable  level 
without  adding  noise  contamination  from  the  measuring 
equipment.  This  difficulty  occurs  because  noise  signals 
are  usually  small  quantities  and  are  sometimes  buried  in 
the  background  noise.  We  designed  and  built  a new  noise 
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Figure  6.8  A plot  of  S-parameter  (S^  ) 


Current  (juA) 


149 


Figure  6.9  I v characteristics  of  the  device  UM516. 
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Figure  6.10  C-V  characteristics  of  the  device  UM516 
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Figure  6.11  A plot  of  S-parameter  (S^ ) 
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measuring  system,  improving  the  sensitivity  over  that  of 
conventional  systems. 

To  measure  a wide  frequency  range  of  the  noise 
spectrum,  we  need  to  employ  several  different  measurement 
schemes.  Basically  we  divided  our  frequency  span  from  100 
Hz  to  1 GHz  into  three  frequency  ranges  because  one 
receiver  can't  cover  the  entire  frequency  span.  In 
addition  it  is  desirable  to  use  different  measuring  set- 
ups to  have  a maximum  sensitivity  for  each  frequency 
range.  All  measurement  set-ups  were  tested  for  linearity 
and  calibrated.  For  calibration,  a commercial  p-i-n 
photodiode  was  used  with  an  identical  device  mount  as  that 
of  the  DUT.  Since  high  frequency  measurements  are  much 
more  difficult  than  low  frequency  measurements  and  require 
state  of  the  art  techniques,  we  will  give  more  attention 
to  the  high  frequency  experiment  than  to  the  low  frequency 
set-ups . 

6.4.1.  High  Frequency  (150  MHz-1  GHz)  Noise  Measurement 
Set-up 

The  highlight  of  this  measuring  system  is  the 
combination  of  a modulation  scheme  first  used  by  Dicke 

[76]  and  the  circulator  method  developed  by  Gasquet  et  al. 

[77]  . Known  as  the  Dicke  radiometer,  the  modulation  scheme 
is  used  to  separate  out  the  weak  signal  buried  in  the 
background  noise  stemming  from  the  thermal  noise  of  the 
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load  resistor  and  the  amplifier  noise.  The  sensitivity  is 
given  by  [78] 


where  B is  the  bandwidth  of  the  receiver  and  T is  the  time 
over  which  the  signal  is  averaged.  In  our  measurements  B 
is  4MHz  and  T has  a value  of  50  seconds.  By  substituting 
these  values  and  using  T = 300  K,  we  obtain  AT  = 0.02K. 
This  is  the  calculated  sensitivity.  In  reality,  available 
sensitivity  is  only  0.2  K,  mainly  because  of  the 
unbalanced  temperature  drift  ratios  between  the 
calibration  source  resistance  and  the  load  resistance.  In 
the  high  frequency  experiment,  impedance  matching  is 
important  because  a false  signal  results  from  an  impedance 
mismatch  during  the  Dicke  switching.  This  impedance 
mismatching  problem  can  be  solved  by  using  a circulator. 
Thus,  regardless  of  the  impedance  changes  in  the  device, 
the  amplifier  is  always  looking  at  a 50-Q  termination 
through  the  circulator.  Several  circulators  are  used  to 
cover  the  wide  frequency  range  (150  MHz  - 1000  MHz)  . 
Detailed  information  about  these  circulator  is  given  in 
Table  6.2. 
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Table  6.2 

Circulators  used  in  the  high-frequency 
noise  measuring  set-up 

Circulator  Name 

Frequency  range 

Loss  (dB) 

Innowave  1M16CV 

147  - 174  MHz 

0.5 

Innowave  1M35CZ 

300  - 400  MHz 

0.4 

P&H  B1-U36310 

500  - 1000  MHz 

0.2 

A block  diagram  of  the  entire  system  is  shown  in  Fig. 
6.12.  An  LED  is  used  as  a light  source  because  of  its  full 
shot-noise  performance.  The  wavelength  of  the  LED  is  1.3 
^im.  To  avoid  stray  light  which  will  cause  parasitic 
absorption,  light  need  to  be  focused  on  the  top  surface 
through  the  window.  Optical  fiber  guiding  turned  out  to  be 
the  best  solution  for  confined  illumination.  The  tip  of  an 
optical  fiber  was  mounted  on  an  XYZ  micrometer  positioner 
and  could  be  easily  positioned  at  the  center  of  the 
window.  A microscope  with  100X  multiplication  is  used  to 
check  the  optical  fiber  alignment  visually. 

The  biasing  circuit  is  shown  Fig.  6.13.  When  we 
increase  the  reverse  bias,  breakdown  will  occur  abruptly. 
Once  breakdown  occurs,  MQW  structures  are  easily 
destroyed.  For  protection  we  built  an  accurate  current 
limiter  whose  resolution  is  as  accurate  as  O.lfiA.  As  a 
power  source,  a regulated  power  supply  is  used.  The  device 
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Figure  6.12 


High-frequency  noise  measurement  set-up. 
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Current  limiting  circuit  for 
the  device  protection. 


Figure  6.13 
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under  test  (DUT)  is  mounted  on  a 50-Q  microstripline  with 
a 50-Q  termination  as  shown  in  Fig.  6.14. 

Immediately  following  the  DUT  is  a coaxial  relay 
switch  which  is  operated  by  a pulse  generator.  The  output 
of  the  relay  switch  is  fed  to  circulator  port  2 while  a 
noise  source  (HP346B)  is  hooked  up  to  circulator  port  1. 
Since  the  HP346B  is  not  a variable  calibration  noise 
source,  an  attenuator  is  used  to  adjust  the  calibration 
noise  signal  level.  With  an  attenuator  in  series, 
calibration  noise  temperature  Tcai  is  given  by  [79] 

TCal=  LTNS  +(1-L)T0  (6.2) 

where  Tns  represents  the  noise  source  temperature,  T0  is 
room  temperature  and  L denotes  the  attenuator  factor. 

Following  the  circulator  is  the  bias  tee  (HP11590A) 
which  couples  the  noise  signal  to  the  low  noise  amplifier 
(LNA)  and  shorts  the  50-Q  load  resistor  D.C.  wise.  The 
latter  helps  to  avoid  Joule  heating  of  the  load  resistor. 

Frequency  selection  is  obtained  with  the  noise  figure 
meter  (HP8970A) , which  is  capable  of  receiving  RF  signals 
from  10  MHz  to  1.5  GHz  with  a bandwidth  of  4.0  MHz. 
Instead  of  using  the  power  detector  built  in  the  noise 
figure  meter,  we  employ  a separate  power  detector 
(HP432B) . In  other  words,  the  noise  figure  meter  is  used 
as  a frequency  converter  in  this  set-up  and  we  get  the 
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Figure  6.14  Configuration  of  a device  mount 
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signal  output  from  its  intermediate  frequency  (21  MHz) 
output.  To  obtain  an  additional  20dB  boost  of  this  signal, 
a high  frequency  amplifier  (HP8447F)  is  placed  after  the 
21  MHz  low-pass  filter.  Power  detector  (HP432B)  produces  a 
voltage  signal  at  the  monitor  termination  as  well  as  an 
analog  reading.  The  voltage  output  is  fed  to  the  FFT 
spectrum  analyzer  (HP3582A)  through  the  low-pass  filter 
with  bandwidth  of  1 Hz.  Finally,  the  FFT  analyzer  displays 
the  magnitude  of  the  detected  power  signal  as  a function 
of  frequency.  When  the  coaxial  relay  switches  slowly 
(0.2Hz)  between  port  1 and  port  2,  the  difference 
between  the  noise  generated  in  the  DUT  and  the  calibration 
noise  source  is  displayed.  The  FFT  spectrum  analyzer 
averages  the  block  function  which  is  collected  for  50 
seconds  (10  periods)  and  separates  out  the  0.2  Hz 
component  by  Fourier  transformation. 

6.4.2.  Medium  High  Frequency  (10  MHz  - 300  MHz)  Noise 
Measurement  Set-up 

In  this  frequency  range,  it  is  difficult  to  obtain 
the  proper  circulators  to  cover  a wide  frequency  range. 
Without  the  circulator,  impedance  mismatches  in  the  front- 
end  cause  a false  signal  during  the  Dicke  switching.  But 
it  turned  out  that  below  300  MHz  the  false  signal  due  to 
the  impedance  mismatching  is  not  serious.  In  that  case,  we 
are  able  to  eliminate  the  offset  caused  by  the  impedance 
mismatches  of  the  front-end  system  by  following  the 
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procedures  which  we  will  discuss  in  the  next  section. 
However,  impedance  changes  due  to  the  bias  dependent 
junction  capacitance  can  not  be  compensated.  Therefore 
this  limitation  allows  us  to  measure  the  noise  of  the  DUT 
at  low  bias  only,  because  then  we  can  neglect  the  change 
in  the  junction  capacitance  with  bias. 

Another  difficulty  arises  when  we  operate  the  relay 
switch  without  using  the  circulator.  For  the  short  time 
interval  when  the  mechanical  switch  moves  from  port  1 to 
port  2,  the  input  of  the  amplifier  floats.  This  results  in 
a large  spike  in  the  output  signal.  To  eliminate  this 
troublesome  spike,  we  designed  a device  whose  block 
diagram  is  shown  in  Fig.  6.15.  The  pulse  generator  which 
operates  the  relay  switch  triggers  the  monostable  multi- 
vibrator (SN74121)  to  generate  the  pulse.  Controlling  the 
potentiometer  we  adjust  the  pulse  width  to  activate  relay 
2 for  a time  interval  which  should  be  slightly  larger  then 
the  duration  of  the  spike.  When  the  relay  is  turned  on, 
the  input  of  the  spectrum  analyzer  floats,  but  the  output 
signal  of  the  spectrum  analyzer  remains  the  same  because 
of  the  large  RC  time  constant  due  to  the  large  value  of 
the  coupling  capacitor  and  the  high  input  impedance.  Thus, 
the  signal  fed  to  the  spectrum  analyzer  is  skipped, 
eliminating  the  spike.  Other  portions  of  the  set  up  are 
the  same  as  in  the  previous  section. 
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Figure  6.15  Measuring  sep-up  with  the  spike-eliminating 
timing  circuit  and  related  waveforms 
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6.4.3.  High  Frequency  Noise  Measurement  Procedure 

The  experimental  procedure  consists  of  four 
measurements  to  determine  the  noise  temperature  of  the 
DUT.  The  following  symbols  will  be  used  in  this  section 
T0  = Room  temperature 

Td0  = Temperature  of  the  DUT  in  equilibrium 
Td  = Temperature  of  the  DUT 
tnd  = Noise  temperature  of  the  DUT 
Ta  = Temperature  of  the  amplifier 
Tcai  = Temperature  of  the  calibration  noise  source 
G2BW  = Gain  bandwidth  product 
r = Voltage  reflection  coefficient 

k = Boltzmann's  constant 

Measurements  proceed  as  follows: 

1)  Short  the  input  and  turn  off  the  calibration  noise 
source.  Then  the  noise  source  temperature  equals  room 
temperature.  This  procedure  gives  a meter  reading  equal  to 

Mi  = G2-BW[kTa  + kT0]  (6.3) 

2)  Keep  the  input  shorted  and  turn  on  the  noise 
source.  Then  the  noise  source  radiates  the  calibration 
noise  power  kTcai  per  unit  bandwidth.  The  meter  reading 
becomes 


M2  = G2-BW[kTa  + kTcai] 


(6.4) 
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3)  Have  the  DUT  connected  to  the  input  and  turn  off 
the  noise  source.  This  gives  a meter  reading  of 

M3  - G2-BW[kTa  + kTDo(l-|r|2)  + kToiri2]  (6.5) 

4)  Have  the  DUT  connected  and  turn  on  the  noise 
source.  The  meter  reading  becomes 

M4  - G2-BW[kTa  + kTood-in2)  + kTcalin2]  (6.6) 

Manipulation  of  the  above  four  equations  yields  the 
following  expression  for  the  DUT  temperature. 

Tdo  = T0[l  + Tcal  (M3  _ Mi)+(M2  - M4)^  (6.7) 

where  M3  - Mi  and  M2  -M4  are  obtained  by  operating  the 
switch  Si  with  the  noise  source  turned  off  and  turned  on, 
respectively. 

Ideally  Tdo  should  be  equal  to  T0.  But  this  is  not 
what  we  experienced.  There  still  might  be  impedance 
mismatches  in  the  front-end  and  small  temperature 
differences  between  the  noise  source  and  the  load  resistor 
in  the  device  mount.  These  non-idealities  result  in  an 
offset  signal  in  the  meter  reading.  We  repeat  the  above 
four  procedures  with  the  DUT  biased.  Then  we  obtain 

(M3  - Mi)  + (M2  - M4)  ] 


Td  — Tq[1  + Tcai 


(6.8) 
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By  subtracting  Eq.  (6.7)  from  Eq.  (6.8)  we  obtain  the  noise 
temperature  Tnd  of  the  DUT. 

Depending  on  the  measuring  frequency  range  we  have 
to  replace  the  circulator. 

6.4.4.  Low-frequency  Noise  Measurement  Set-up  and 
Procedure 

As  shown  in  Fig.  6.16,  in  this  measurement  set-up  a 
standard  measurement  scheme  is  used.  Sensitivity  is  not  a 
critical  limitation  in  this  frequency  range  because  of  two 
reasons.  One  is  that  employing  a large  value  (5  KQ)  load 
resistor  improves  the  sensitivity  hundred  times  over  the 
high-frequency  set-up  value.  The  second  reason  is  that 
preliminary  experiments  showed  large  excess  noise  levels. 

Since  the  frequency  range  of  the  LNA  and  spectrum 
analyzer  is  not  wide  enough  to  cover  the  entire  frequency 
range  from  100  Hz  - 1.5  MHz,  required  in  this  system,  we 
divide  this  frequency  range  into  two  frequency  spans 
ranging  from  100  Hz  - 25  KHz  and  from  10  KHz  - 1.5  MHz. 
For  each  frequency  span,  the  required  equipment  is  listed 
in  Table  6.3. 


Table  6.3 


100  Hz  - 25  KHz 

10  KHz  - 1.5  MHz 

Analyzer 

HP 

Wave  Analyzer 

LNA 

EG&G 

ANZAC  or  CHW 

Noise  Source 

Signal  supplied 
by  HP  analyzer 

Rhode  & Schwarz 
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Figure  6.16  Low-frequency  noise  measuring  set-up 
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Three  measurements  are  required  to  determine  the  noise  of 
the  DUT.  The  procedure  is  as  follows. 

1)  Short  the  input  of  the  LNA.  Then  the  meter  reading 
corresponds  to  (see  Fig.  6.17a) 

Mi  = G2BW-SVa  (6.9) 

where  Sva  represents  the  voltage  spectral  density  of  the 
LNA.  Note  that  squared  values  of  reading  Mi  and  gain  G are 
used  since  noise  is  a quadratic  quantity. 

2)  Connect  the  DUT  to  the  LNA.  The  meter  reading 
becomes  (see  Fig.  6.17b) 

M^  - G2BW-[SVa  + (SiL  + SiD  + Sia)  (Ri//RL)  2 (6.10) 

where  Sil  = current  spectral  density  of  the  load  resistor 

Sio  = current  spectral  density  of  the  DUT 

Sia  = current  spectral  density  of  the  LNA 

Ri  = internal  impedance  of  the  LNA 
Rl  = load  resistance 

3)  Add  the  calibration  noise  source  with  the  DUT 
connected  to  the  LNA.  This  gives  a meter  reading  of 

M^  = G2BW-  [SVa  + (SiL  + SiD  + Sia  + SCal)  (Ri//RL)2] 

(6.11) 

Note  that  Fig.  6.17c  shows  the  equivalent  circuit 
concerning  the  above  equation. 


167 


Figure  6.17  Noise  equivalent  circuit  of  the  measuring 
set-up 
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Manipulating  these  four  equations  yields  the  unknown 
spectral  density  of  the  DUT 


SiD 


- M? 
M3  “ 


Sycal 


Rt 


4KT 

Rl 


>ia 


(6.12) 


There  is  an  alternative  scheme  which  simplifies  the 
previous  procedure.  Instead  of  using  a calibration  noise 
source,  the  thermal  noise  of  the  load  resistor  can  be  used 
to  calibrate  the  measurements.  In  the  frequency  span  from 
10  KHz  - 1.5  MHz,  this  alternative  scheme  was  used.  Both 
of  the  schemes  yielded  the  same  results.  Following  is  the 
modified  procedure. 

1)  Short  the  input  of  the  LNA 

2)  Connect  the  DUT  with  bias  turned  on 

3)  Keep  the  DUT  connected  and  turn  off  the  bias 


These  procedures  yield  the  following  three  equations 

M^  = G2BW- [SVa]  (6.13) 

M^  = G^BW- [SVa  + Sth,RL  + SiD(RL//Ri)2]  (6.14) 
M3  = G2BW- [SVa  + Sth,RLl  (6.15) 

From  these  three  equations,  we  obtain 

m|  - m|  4KT 

iD  = 


(6.16) 
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6.5.  Results  and  Discussion  of  the  Noise  Measurements 
The  noise  measurements  on  device  UM505  were  carried 
out  first.  The  device  was  kept  at  room  temperature.  The 
results  for  the  current  noise  spectra  between  100  Hz  and  1 
GHz  are  presented  in  Fig.  6.18  We  note  that  the  DUT 
exhibits  excess  noise  which  decreases  approximately  as  1/f 
at  low  frequencies  and  reaches  a plateau  level  for  high 
frequencies.  For  an  ordinary  photodiode  the  plateau  level 
is  supposed  to  coincide  with  full-shot  noise.  However, 
this  MQW-APD  exhibits  sub-shot  noise  at  low  bias.  For 
example,  for  a bias  current  of  4 |1A,  the  measured  current 
noise  is  a factor  of  four  lower  than  the  full-shot  noise 
level.  As  the  bias  increases,  the  noise  increases  and 
full-shot  noise  is  obtained  at  a bias  current  of  9 HA.  To 
investigate  the  bias  dependence  of  the  noise  in  more 
detail,  several  measurements  were  carried  out  at  different 
frequencies.  Their  results  are  presented  in  Fig  6.19  to 
Fig.  6.22.  In  these  measurements,  the  bias  range  is 
between  3 |iA  and  10  |1A.  The  lower  limit  was  determined  by 
the  maximum  sensitivity  of  the  measuring  system  and  the 
upper  limit  is  chosen  to  avoid  excess  tunneling  current. 
To  find  the  noise  level  at  low  bias  below  the  measuring 
range,  we  used  an  extrapolation  method  and  obtained  0.06  K 
for  1.0  }1A,  which  is  approximately  1/60  smaller  than  full- 
shot  noise.  This  sub-shot  noise  phenomenon  is  a result  of 
the  trapping  and  detrapping  in  the  undepleted 
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Noise  Spectrum  (Device:  #UM505;  90  — 90  A) 


Figure  6.18  Noise  spectra  of  the  device  UM505. 
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Figure  6.19  Noise  temperature  vs.  bias  current 

at  an  operating  frequency  of  150  MHz 
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Figure  6.20  Noise  temperature  vs.  bias  current 

at  an  operating  frequency  of  300  MHz 
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Figure  6.21  Noise  temperature  vs.  bias  current 

at  an  operating  frequency  of  400  MHz 


Noise  Equivalent  Temperature  (K) 


174 


1 


10 


UM505 

Freq  = 700  MHz 


0 


0 


i i i i i i 


10°  Bias  Current  (|1A) 


Figure  6.22  Noise  temperature  vs.  bias  current 

at  an  operating  frequency  of  700  MHz 
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multiquantum-wells  as  mentioned  in  chapter  V.  From  the 
noise  spectrum,  we  also  note  that  the  large  excess  noise 
levels  showing  up  over  a wide  frequency  range  need  to  be 
explained.  One  possible  explanation  is  that  the  large 
number  of  heterostructure  interfaces  results  in  a large 
number  of  interface  states  or  deep  level  traps  which  will 
increase  the  generation-recombination  process. 

The  noise  spectrum  of  device  UM516  was  measured 
between  10  Hz  and  25  KHz.  Fig.  6.23  demonstrates  the 
results  and  we  see  that  this  device  also  exhibits  1/f- 
noise  like  excess  noise  at  low  frequencies.  Fig.  6.24 
presents  the  high  frequency  noise  spectra  which  were 
measured  between  10  MHz  and  1 GHz  at  room  temperature.  The 
noise  levels  off  at  roughly  150  MHz.  Below  this  corner 
frequency,  the  frequency  dependence  is  weak  in  contrast  to 
what  we  observed  in  device  UM505.  Our  explanation  is  as 
follows.  Device  UM516  has  25  periods  of  heterostructure 
layers  compared  to  100  layers  for  UM505.  Having  a smaller 
number  of  interfaces  results  in  a smaller  number  of 
interface  states  or  traps,  thus  resulting  in  a decrease  in 
generation-recombination  noise. 

The  current  noise  versus  bias  current  was  measured  at 
several  different  frequencies.  The  bias  range  is  between  4 
|IA  and  20  jiA.  Above  7 JiA,  the  noise  of  the  DUT  is 
increasing  nearly  linearly  with  the  bias  current,  but  is  a 
factor  of  three  less  than  the  full-shot  noise  level.  Even 
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Figure  6.23  Low-frequency  noise  spectra 
of  the  device  UM516 
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Figure  6.24  Noise  spectra  of  the  device  UM  516. 
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at  high  bias,  this  device  exhibits  sub-shot  noise.  The  C-V 
characteristic  also  indicates  that  punch  through  or  full 
depletion  is  difficult  to  obtain  in  this  structure.  This 
suggests  that  there  are  at  least  three  or  four  undepleted 
quantum-wells  in  the  center  of  the  depletion  region. 
Compared  to  device  UM505,  this  structure  (UM516)  has 
wider  quantum-wells  (500  K)  and  a larger  number  of 
thermally  generated  carriers.  When  we  increase  the  bias 
further  to  deplete  the  undepleted  quantum-wells,  large 
tunneling  currents  in  the  low  band-gap  may  prevent 
complete  depletion. 

In  Fig.  6.25  a break  point  occurs  at  bias  current  of 
7 |1A.  Below  this  bias,  it  is  difficult  to  estimate  the 
factor  of  sub-shot  noise  due  to  the  poor  sensitivity  of 
the  set-up.  To  overcome  this  difficulty,  we  used  a 
semiconductor  laser  as  a light  source  to  get  more  photo- 
current. The  noise  of  the  laser  was  characterized  with  a 
conventional  p-i-n  diode.  Fig.  6.26  presents  the  bias 
dependent  noise  of  the  DUT  illuminated  by  the  laser.  At  15 
MA  bias,  which  is  the  lower  bias  limit,  the  device  noise 
is  1/23  smaller  than  full-shot  noise.  The  reciprocal  value 
of  this  fraction  is  close  to  the  number  of  quantum-wells 
in  this  structure. 

Device  UM504  has  very  thin  (30  A) , but  a large  number 
(180)  of  heterostructure  layers.  Therefore,  low-frequency 
noise  measurements  exhibit  large  excess  noise  level  as 
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Figure  6.25  Noise  temperature  vs.  bias  current 

at  an  operating  frequency  of  500  MHz 
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Figure  6.26  Noise  temperature  vs.  bias  current 

at  an  operating  frequency  of  700  MHz 
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shown  in  Fig.  6.27.  The  high  frequency  noise  spectrum 
shown  in  Fig.  6.28  seems  to  indicate  that  the  current 
noise  of  the  DUT  is  strongly  frequency  dependent.  But  if 
we  take  the  series  resistance  and  junction  capacitance 
into  account  and  calculate  the  actual  amount  of  noise 
generated  in  the  DUT,  the  spectrum  becomes  flat  and 
coincides  with  a full-shot  noise  level.  This  has  to  be 
expected  because  the  structure  with  thin  transparent 
layers  is  quite  similar  to  a conventional  bulk  structure. 
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Figure  6.27 


Noise  Spectra  of  the  device  UM504 
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Figure  6.28  High-frequency  Noise  Spectra 
of  the  Device  UM504 


CHAPTER  VII 

SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FUTURE  RESEARCH 

7 . 1 Summary  and  Conclusions 

We  have  described  a theoretical  and  experimental 
study  on  the  opto-elect ronic  properties  of  MQW 
heterostructure  p-i-n  diodes  grown  by  molecular  beam 
epitaxy. 

In  the  course  of  our  experimental  investigations,  we 
have  faced  many  unexpected  difficulties.  First  of  all,  the 
large  excess  noise  components  at  low  frequencies  resulted 
in  a lower  corner  frequency  of  the  high  frequency  noise 
plateau  of  200  MHz,  which  is  an  extraordinarily  high  value 
compared  to  the  value  of  several  MHz  in  conventional  p-i-n 
diodes.  Hence,  we  had  to  build  a new  high-frequency  noise 
measuring  set-up  employing  Gasquet ' s circulator  method  and 
Dicke's  radiometer  technique. 

Secondly,  in  our  high  frequency  measurement,  the 
noise  signal  was  often  masked  by  parasitic  effects  arising 
from  junction  capacitance,  series  resistance,  and  signal 
cables.  Therefore,  S-parameter  measurements  were  performed 
to  characterize  the  high-frequency  small  signal  properties 
of  the  device  and  the  device  holder,  so  that  we  could 
correct  the  measured  noise  signal  for  parasitic  effects. 
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Thirdly,  distinctive  sub-shot  noise  levels  observed 
in  the  devices  with  90-90  and  500-500  Angstrom  layers  were 
not  understood  at  first.  To  confirm  the  credibility  of  the 
measuring  system,  the  noise  of  a conventional  p-i-n  diode 
mounted  in  an  identical  device  mount  was  measured.  The 
results  of  this  measurement  assured  us  that  our  measuring 
system  was  working  properly. 

As  suggested  by  Gisolf  [65],  Bernamont  [66],  and 
Machlup  [67] , a possible  cause  for  distinctive  sub-shot 
noise  in  semiconductor  devices  can  be  the  coupling  of 
carrier  motion  to  the  external  circuit  when  the  carriers 
are  trapped  and  detrapped.  We  found  a strong  correlation 
between  the  number  of  undepleted  quantum-wells  and  the 
sub-shot  noise  levels  over  a wide  bias  range.  By 
incorporating  the  results  of  C-V,  I-V  and  noise 
measurements  as  a function  of  bias,  we  conclude  that  the 
observed  sub-shot  noise  levels  are  due  to  the  trapping  and 
detrapping  in  the  undepleted  multiquantum-wells.  The 
following  observations  can  be  understood  in  this  context. 

In  the  case  of  the  30-30  Angstrom  devices,  which  have 
transparent  quantum-well  layers,  full-shot  noise  was 
observed  in  the  frequency  range  between  200  MHz  and  1 GHz, 
whereas  in  the  case  of  the  90-90  and  the  500-500  Angstrom 
devices,  where  the  barriers  are  no  longer  transparent  and 
the  carriers  are  trapped  and  thermal  assisted  detrapped  in 
the  quantum-wells,  distinctive  sub-shot  noise  levels  were 
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observed.  As  reverse  bias  increases,  the  quantum-well 
layers  located  in  the  i-region  are  depleted  gradually 
i.e.,  one  quantum-well  at  a time,  resulting  in  a gradual 
increase  of  the  noise  level.  At  the  field  strength  where 
all  quantum-wells  are  depleted,  full-shot  noise  is 
observed. 

The  excess  noise  observed  at  frequencies  below  200 
MHz  is  of  the  generation-recombination  type  and  is  caused 
by  electrons  trapping  and  detrapping  in  heterostructure 
interface  states. 

To  support  the  experimental  study  of  carrier  trapping 
and  detrapping  noise,  a qualitative  theory  was  developed 
based  on  a simplified  MQW  model.  Only  two  energy  levels 
were  assumed  in  this  model,  namely  the  extended  state 
level  which  corresponds  to  the  conduction  band-edge  of  the 
wide-gap  InAlAs  layers  and  the  localized  state  level  which 
corresponds  to  the  conduction  band-edge  of  the  InGaAs 
layers.  To  obtain  the  current  fluctuation  expression,  the 
transfer  impedance  method  was  used  in  solving  the 
transport  equations  to  which  we  added  rate  fluctuation 
noise  terms. 

From  C-V  measurements,  successive  depletion  of 
quantum-wells  over  a wide  bias  range  was  observed  in  90-90 
and  500-500  Angstrom  devices,  whereas  in  30-30  Angstrom 
case,  immediate  depletion  near  zero  bias  was  observed. 
This  C-V  data  in  conjunction  with  bias  dependent  noise 
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data  strongly  indicate  that  trapping  and  detrapping  in  the 
undepleted  quantum-wells  determine  the  overall  current 
noise  level. 

The  dark  current-voltage  characteristic  is  controlled 
by  field  assisted  thermal  emission  of  thermally  generated 
carriers  from  the  quantum-wells,  whereas  the  photocurrent 
is  controlled  by  field  assisted  injection  of  optically 
generated  electrons  from  the  low  bandgap  absorption  region 
into  the  depletion  region  at  low  bias.  At  higher  bias,  the 
photocurrent  is  limited  by  carrier  transit  across  the 
undepleted  quantum-well  layers.  At  extremely  high  bias, 
Zener  tunneling  current  becomes  dominant. 

7.2.  Recommendations  for  Future  Research 

Firstly,  we  recommend  more  detailed  C-V  measurements 
to  investigate  the  charge  distribution  of  the  MQW 
structures.  From  these  measurements  accurate,  bias 
dependent,  band  profiles  can  be  obtained. 

Secondly,  we  recommend  that  structures  made  of  other 
compound  materials  and  different  compositions  should  be 
investigated  to  compare  the  characteristics  of  charge 
transport  through  the  MQW  layers. 

Thirdly,  bottom  illumination  is  recommended  to  give 
pure  hole  injection  instead  of  electron  injection  as  was 
employed  in  our  study.  In  this  way,  hole  transport  through 
the  multiquantum-well  hetero- junctions  can  be  compared 
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with  electron  transport  in  terms  of  detrapping  activation 
energy,  bias  dependence  of  hole  injection  over  a valence- 
band  potential  barrier,  hole  confinement  in  quantum-wells, 
etc. 

Fourthly,  we  recommend  employing  a filter  with  a 
wider  bandwidth  in  the  high-frequency  receiver  to  enhance 
the  measuring  sensitivity.  Currently  our  receiver  has  a 4 
MHz  bandpass  filter.  For  instance,  with  the  use  of  16  MHz 
bandpass  filter,  sensitivity  will  be  enhanced  two  times, 
allowing  us  to  get  more  data  in  the  lower  bias  regime. 

Fifthly,  low  temperature  noise  measurements  are 
recommended  to  accurately  determine  the  activation  energy 
of  the  detrapping  carriers  in  the  quantum-well  layers. 
Measurements  of  the  low-frequency  spectrum  as  a function 
of  temperature  will  greatly  benefit  the  analysis  and  the 
identification  of  an  excess  noise  source.  Also  the  dark 
current  can  be  appreciably  reduced  by  lowering  the 
temperature.  Then,  it  might  be  possible  to  increase  the 
reverse  bias  to  initiate  band-to-band  ionizations.  One  of 
the  practical  problems  that  needs  to  be  addressed  is  how 
to  put  an  optical  alignment  mechanism  inside  the  cryostat 
chamber. 

Sixthly,  we  recommend  that  for  low-noise  APD 
applications,  the  desired  enhancement  of  ionization 
coefficient  ratio  a/p  should  be  obtained  by  employing 
resonant  impact-ionization  schemes.  Hence,  HgCdTe  or 
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GaAlSb  bulk  materials  with  split-off  bands  seem  to  be 
relevant.  From  our  study  on  trapping  and  detrapping  noise, 
it  is  found  that  homogeneous  bulk  structures  are  to  be 
preferred  since  they  have  less  excess  noise  at  low- 
frequencies  . 

Finally,  we  recommend  that  Monte  Carlo  simulations 
should  be  performed  to  gain  more  understanding  of  the 
carrier  transport  mechanism  in  the  direction  perpendicular 
to  the  hetero junction  interface. 


APPENDIX  A 

APPLICATION  OF  THE  COMPOUNDING  THEOREM 
TO  A BERNOUILLI  TRIAL 

Consider  the  following  situation: 


Bemouilli  Trial 


/ / / / 


The  generating  function  of  the  Fowler  type  of  y and  0y  is 
by  definition  equal  to 

<J)y  (z)  = X zYW  (y) , (A.  1) 

y 

where  W(y)  is  the  probability  that  the  statistical 
variable  y has  a value  y. 

Hence  in  the  above-pictured  situation, 


(t)y(z)  = X zYW(y)  = z°W(y=0)  + z1W(y=l)  + z2W(y=2)  + • • • 

y 

(A. 2) 

The  only  way  that  y can  be  zero  is  by  not  feeding  the 
statistical  process  a particle.  The  probability  for  this 
is  1 - X;  consequently,  W(y=0)  = 1 - X.  From  inspection  of 
the  above  figure,  we  conclude  that  W(y=l)  has  to  be  equal 
to  A.-W(u=l),  and  that  W(y=2)  equals  A,*W(u=2),  etc. 
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Substituting  these  results  in  equation  (A. 2)  gives  for 

<|)y  ( z ) : 

<])y  (z ) = l-(l-X)  + ziX,W(u=l)  + z2A,W(u=2)  + ••• 

= 1 - X + ^(|)u(z)  (A. 3) 

If  we  define  the  function  co(z)  by 

co(z)=  1 - X + A,z, 

then,  Eq.  (A. 3)  reads 

(j)y(Z)  = co(<J)u  ( z ) ) . (A.  4) 

Equation  (A. 4)  is  a form  of  the  compounding  theorem  we 
used  in  Chapter  II,  Eq.  (2.94). 


APPENDIX  B 

EVALUATION  OF  SUB-INTEGRALS 


xi 

/»  x ' 

Jpdx" 

Q (x' ) e°  dx' 


J 

o 


Qi  [e-Piii- 


pi 


i] 


x2 

/»  x ' 

Jpdx" 

Q(x' )e0  dx' 


J 

XI 


XI  X ' 

Jpdx"+  Jpdx" 
eo  xi  dx' 


J 

XI 


X2 

Q2  J -ep2  ^x  _Xl^dx' 

xi 

X2 

Q2 -e^1^1 -e~^2Xl  J ep2x  dx' 

xi 


Q2.e(pi"p2)1i.^-[ep2x2  - e"p2xi] 


|2  .e(P1-p2)li[eP2(li+l2)  _ g^1! 


QZ.  j-gP ili+P 2I2  _ epili] 


193 
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